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ABSTRACT 


This paper presents a method for estimating the forces and moments on the six 
axes of an Unmanned Underwater Vehicle (UUV) operating in close proximity to a host 
vessel. Forces and moments encountered by the vehicle in the host vessel launchway, near 
the vessel and clear of the launchway, and far away from the host vessel, are modeled 
from experimental and theoretical sources. A quasi - steady approach to the vehicle 
launch dynamics allows coupling of experimental data, potential flow theory, and 
hydrodynamic maneuvering coefficients to develop an algorithm capable of modeling 
vehicle and host vessel interactions with any user defined trajectory. 


A vehicle motion simulation study was performed using the Revised Standard 
Equations of Motion (EOM) augmented to incorporate the interaction force model for 
proximity effects of the UUV to a host. Application of the trajectory simulation and force 
and moment time history are developed for a underwater vehicle launching from a 
submerged host illustrating the hydrodynamic interaction of the vehicle to the host. The 
method is developed to account for host forward speed, UUV propulsor effects, and open 
water hydrodynamics. The robustness of the method is demonstrated for several launch 
trajectories. 
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Chapter I INTRODUCTION 


1.1 Motivation 


Underwater vehicles find applications in an increasing number of diverse missions 
and roles. Recent uses of AUV's and ROV's both in commercial and military 
applications continue to outline the capabilities, and limitations of employing unmanned 
underwater vehicles (UUV). Integral to the use of a UUV, whether autonomous 
underwater vehicle (AUV) or remotely operated vehicle (ROV), is the process used for 
launch and recovery. This process is largely driven by the type of support ship, land based 
launch and recovery facility, or submerged host platform. Launch and recovery from a 
submerged host represents an opportunity to add flexibility to the vehicle mission, but 
adds complexity to vehicle design. The ability of a submerged host vessel, hereafter 
generically referred to as the host, to deploy a ROV or AUV to a mission commencement 
location in latitude, longitude and depth offers unique advantages in mission duration, real 
time interrogation, and covert mission accomplishment over other launch platforms. The 
ability to properly control a vehicle as it transitions from within the submarine launchway 
to a point far away from the host and back until finally rehoused represents a complex 
dynamic process for which improved models are necessary to investigate the forces, 
moments and trajectory of the vehicle during the course of its mission and through launch 
and recovery. The design of effective vehicle control is the goal of the hydrodynamic 


modeling of the launch and recovery process. 








1.2 Research Objectives 


Launch and recovery of a vehicle from a submarine is expected to represent the 
most demanding case for control effectiveness. The effect of the host body speed and 
launchway geometry on the hydrodynamic flow characteristics in and around the 
launchway have been documented by experiment [11]. Additionally, interaction effects 
between the UUV and the host have been quantified by experiment for specific cases [11]. 
One can foresee that the combined effects which in turn are coupled, lead to a complex 
hydrodynamic modeling problem. 

The first objective of this thesis is to use experimental observations from a scale 
UUV within and in the vicinity of a launchway to develop a model to describe the forces 
and moments acting on a UUV as a function of hydrodynamic and geometric criteria 
relative to the host. Due to the specific nature of vehicle and host geometry on the forces 
experienced, the model will be tailored to the host and UUV. However this specificity is 
necessary only in so far as to proceed to the next objective, that of determining a 
reasonable simulation to determine the time dependent forces and moments acting upon 
the vehicle. 

This second objective is to incorporate this force model into a form of the Revised 
Standard Submarine Equations of Motion [5] to determine the forces and moments acting 
upon the vehicle given a specified trajectory. This investigation relies upon mating 
experimental data in the complex region within the launchway, with inviscid flow code 
results external to the launchway and in proximity to the host "boundary", and accounting 
for viscous effects of cross flow drag associated with the vehicle transitioning through a 
region where the vehicle axis is not aligned with the mean free stream velocity. The 
resulting control force simulator is intended to be adaptable to alternate vehicle 


geometries. 








Chapter2 FORCE MODEL 


2.1 Introduction 


The forces acting upon a UUV operating in the open ocean differ from those 
experienced when operating in the vicinity of a boundary, such as the free surface, ocean 
bottom, or host vessel. For the purposes of this thesis, boundary interaction effects are 
limited to a plane rigid wall approximation of a host vessel except in the internal region of 
the launchway ofthe host. Interaction forces during launch and recovery, and through the 
transition to the "infinite fluid", are of direct interest in determining the feasibility of the 
control system and/or the operating envelope for launch and recovery. Significant work 
has been undertaken in the past, including computational modeling of these forces to 
assess the effects of high speed ejection of torpedoes from submarines [2]. In the region 
of interest for UUV launch and recovery, it will be necessary to operate within a different 
regime of slow forward speed of the host and slow relative speed of the UUV to the host. 
This combination of slow host and UUV relative speed indicates the plausibility of 
modeling the UUV launch and recovery interaction forces as quasi - steady, the UUV 
stationary with respect to the host and slowly changing with respect to its variation in 
geometric relation to the host. This approach allows use of experimental results to 
develop a model for the forces in this region. Steady flow potential based codes may then 
be used to model the interaction effects of the vehicle near the host to resolve areas where 


experimental data is incomplete. 


2.1.1 Background 
The determination of the hydrodynamic forces present on a vehicle transitioning 
from a submarine launchway is complicated by the number of parameters involved and the 


complexity of the fluid -body interactions. For a baseline case with no vehicle in the 
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launchway Sjoblom and Schwemin [15] showed that the flow in a typical submarine 
launchway is unsteady and turbulent, with turbulence intensities within the cavity ranging 
from 10 - 20% of free stream velocity. Experimental observations document average flow 
magnitudes 0 - 1096 of the free stream at the aft end ofthe cavity and approximately 3096 
of free stream at the midpoint of the launchway shutter. They document a large 
recirculation region in mid - region of the cavity and a separated region approximately 
3096 of the cavity length forward. As reported, the turbulence shows a flow periodicity 
with frequency content between 1 and 4 Hz. External to the launchway and turbulent 
boundary layer potential flow predictions compare well with measurements. 

In order to assess the force and moment characteristics of a vehicle in launch or 
recovery transition, both the region within the launchway and external to it require 
modeling. The model must be sophisticated enough to capture three dimensional 
geometry effects, body - to - body interactions, as well as viscous flow effects. Inviscid 
potential flow models capture the proximity interaction effects; however, they do not 
include viscous drag effects including cross flow drag or lift forces generated by control 
surfaces and vorticity shed into the wake from separation. Numerical potential and 
viscous flow techniques can not at present handle all these effects so a combination of 
experimental results, numerical techniques, and semi-empirical non linear motion 


derivatives are suggested for investigating the launch and recovery envelope. 


2.1.2 Experimental Data 


Launchway Zone 

Jensen and Milgram [11] describe an experiment conducted using a 1/10 scale 
model of an UUV transitioning through a submarine launchway. Figure 2-1 illustrates the 
basic configuration and geometry associated with the vehicle and the host. Principle 


dimensions of the UUV vehicle are listed in table 2-1. 
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| UUV Principle Dimensions 


Length (ft) 


diameter (ft 





Table 2-1 Principle Vehicle Dimensions 
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Figure 2-1 


Data was obtained for the forces and moments at various points of advance along the 
launchway centerline. Reference to a "guide can" structure existing in the torpedo tube 
launchway is made in the model tests. The vehicle was advanced from a position farthest 
aft, corresponding to the aft end of the propulsor duct flush with the forward edge of the 
guide can, to a position 10 inches forward of this reference (100 inches full scale). This 
range of traverse starts with the vehicle fully within the launchway and extends to a 
position with the aft end of the propulsor duct approximately even with the forward most 
point of the launchway. Data in all six degrees of freedom were obtained at each position 
along the launchway for different onset flow conditions corresponding to the speed of the 
free stream past the host allowing the effect of forward speed on the host to be modeled. 
The UUV model was configured with a ducted propulsor which was operated during the 
tests at various RPM for each advance position and onset flow combination to simulate 


the influence of propulsor forces and race effects on the induced flow. 
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Near Boundary Zone 

Once the vehicle is outside of the launchway, the geometry of the vehicle and the 
host can be approximated as a body in proximity to a plane rigid boundary. The approach 
taken here, utilizes results from a Boundary Integral Element Method (BIEM) panel code 
to determine the force and moment effects due to proximity operations external to the 
launchway. This method uses inviscid potential flow theory and does not include viscous 
effects. A description ofthe method and implementation is contained in [14]. The effect 
of the propulsor on the induced flow, and accompanying modification to the pressure 
integral is not captured by the BIEM in this application. An attempt to include these 
effects is made through use of experimental data obtained for the UUV at various onset 


flow speeds, propeller RPM, and clearance from the approximated plane boundary. 


2.2 Modeling Results 


2.2.1 Coordinate System 

The force and moment models are referred to a right-hand orthogonal coordinate 
system of moving axes, fixed in the UUV with the origin located at the center of mass 
(CG) of the vehicle. The positive directions of the axes are x - forward, y - starboard, and 


z - downward in accordance with Standard Nomenclature [16]. 


2.2.2 Launchway Model 
The experimental data lead to a model of the forces and moments experienced by 
the vehicle as a function of specific dependent variables. The parameters varied during the 
testing used to develop the model were: 
(1) onset flow speed 
(2) advance along the launchway axis 


(3) propulsor rpm 





(4) thruster flow and direction 

The data runs with the thruster running were not used in this analysis in order to limit the 
force and moment effects to those induced by body and launchway geometry and 
propulsor effects. The forces and moments may be developed in accordance with the 
following dimensional form: 

F = (a uU py) 

= (а.о у) 
The following non-dimensional terms are defined for use in developing a model 
independent of scale. 


Advance 
Length 


CFY - — 

2 

a 
TS 

1 2 

= ( PROPULSIVE FORCE ) 


CMZ = 


РГЕ 7 


where: CF = Force Coefficient 
CM = Moment Coefficient 
PTF = Propulsive Thrust Factor 
Uo, = Onset flow speed 
A, = Frontal Area 
L = Vehicle Length 


For each data set the forces and moments recorded were converted to standard body 
coordinates with forces and moments limited to the horizontal plane being modeled in this 


case. A propulsion model was required which would allow robust modeling of propeller 
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effects. The non-dimensional Propulsive Thrust Factor (PTF) was used, with PTF 
determined for each data set by measuring the mean recorded thrust at zero RPM for each 
onset flow condition and advance to length position and subtracting this value from the 
axial force. The zero PTF data corresponds to the axial drag force experienced in the 
launchway due to the body position alone and was modeled as well. The procedure 
selected to arrive at a model recognized that significant variance existed in the 
experimental data. Figures 2.2 through 2.14 reflect the non dimensional data with 


associated best least squares linear regression fit. 


Drag Force CFX vs A/L [(xposH)A] 


CFX 


eqn: y=a+b*(x-1.47), R:0.107, max dev:0.375 
a=-0.637341, b=-0.853415 





АЛ, 


Figure 2.2 Drag Force 
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СЕҮ(РТЕ) 


СЕҮ(РТЕ) 


0.5 


СЕҮ(РТЕ) (ог АЛ. 0.68 


y =+0.0322x" +0.919, max dev:0.216, r”=0.200 





Figure 2.3 Lateral Force for A/L = 0.68 
СЕҮ(РТЕ) (ог АЛ. 0.84 


y =+0.0773x" +1.05, max dev:0.336, r”=0.210 





РТЕ 


Figure 2.4 Lateral Force A/L = 0.84 
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СЕҮ(РТЕ) 


СЕҮ(РТЕ) 


СЕҮ(РТЕ) (ог АЛ. 1.0 





Figure 2.5 Lateral Force A/L = 1.0 
CFY(PTF) for АЛ, 1.16 





Figure 2.6 Lateral Force A/L = 1.16 
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СЕҮ(РТЕ) 


СЕҮ(РТЕ) 


СКУ(РТЕ) Гог АЛ. 1.32 





Figure 2.7 Lateral Force A/L = 1.32 
СЕҮ(РТЕ) (ог АЛ. 1.48 





Figure 2.8 Lateral Force A/L = 1.48 
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CMZ(PTF) 


CMZ(PTF) 


СМ2ЦРТЕ) (ог АЛ. 0.68 


--0.113х" 40.0339, тах 4еу:0.0745, г-0.971 





Figure 2.9 Lateral Moment A/L = 0.68 
CMZ(PTF) for A/L 0.84 


M -0.0300x' +0.0401, max dev:0.171, r”=0.184 





Figure 2.10 Lateral Moment A/L = 0.84 
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CMZ(PTF) 


CMZ(PTF) 


CMZ(PTF) for A/L 1.0 


y =+0.0656x" +2.37E-4, max dev:0.551, r”=0.397 





Figure 2.11 Lateral Moment A/L - 1.0 
CMZ(PTF) for A/L 1.16 





Figure 2.12 Lateral Moment A/L = 1.16 
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CMZ(PTF) 


CMZ(PTF) 


CMZ(PTF) for A/L 1.32 





Figure 2.13 Lateral Moment A/L = 1.32 
CMZXPTF) for A/L 1.48 


y =+0.191x" -0.0505, max dev:0.0651, r”=0.988 





Figure 2.14 Lateral Moment A/L = 1.48 
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А Опе Way Analysis of Variance (ANOVA) was conducted on each data set of force and 


moment as a function of A/L and PTF. Table 2-2 summarizes the results. 





AL | Data | FRatio | Prob>F | r2. | 
| 
le c | m CM Zona 18539: 60046 ks: 0,000 e 
mmn | m CMZI IEA 
2 | СМ2 | 447450 | «0.0001 | 0.396869 - 
k | 
| 
(ағын 
1.47 













0.184186 


Oia a 


116.0631 | «0.0001 | 0.773429 | 
[0.757141 | 


0.988475 | 


а 






Ex. 
EX. 
CMZ 


C 










| 31.1761 
857.6477 











«0.0001 
T (2) outlying points omitted 
Table 2-2 One Way ANOVA Results 
In Table 2-2 , the F-ratio is a statistical measure indicating whether the model is a 


better fit than the response mean, as determined by 


(У Ds | Degrees of freedom) | 


| 2x ур / Degrees of freedom) | 


The r? category represents the amount of total variation accounted for by the model with 


F -ratio— 


the remaining variation attributed to random error and is determined as 


5 (Sy?) 


r? = model 


2 уг Ctotal 


To conclude from Table 2-2 that the force and moment are dependent upon PTF the F- 
ratio must exceed an observed significance probability (P-Value) of 0.05. To determine 
the effectiveness of the model, r? is used in conjunction with the F-ratio. The case for 

A/L of 1.0 utilized the r parameter to omit two outlying points from the model. The 


results indicate that for the majority of sample conditions the effect of PTF at each A/L is 
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present, and that a linear model is justified. Higher order models were not justified In 
order to arrive at a smooth model for force and moment as a function of PTF and A/L the 
éCFY(PTF)]cPTF ,,, is used to fit a 3rd order cubic function and the intercept is adjusted to 
make a best fit. Figure 2-15 and 2-16 reflect the best fit for force and moment within the 
launchway. It can be seen that a positive outboard force exists for a minimum A/L of 0.64 
and remains positive throughout the vehicle transition along the launchway centerline with 
a maximum value at an A/L of 1.15 which corresponds to approximately one quarter of 
the stern shadowed by the launchway. The lateral force then diminishes to 50% of this 
force at the maximum A/L of 1.48. Lateral moment effects are quite different, showing an 
inboard moment at minimum A/L followed by a zero crossing at A/L — 0.9. From this 
position the moment shows a continual increase to a maximum with the vehicle fully 
emerged from the host, corresponding in magnitude but opposite sign its value at 
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ӘСМ?АРТЕУӘРТЕ 
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y=-2.04x" +5.59x" -4.49x" +1.13, R:0.0260, max dev:0.0212 


0.75 1.00 1.25 


АЛ, 


Figure 2.15 Fit of AFY(PIF)/@PTF , , 
ОСМ2ЦРТЕ)УФРТЕ УБ АЛ, 


y=-0.153x" +0.198x" +0.516x" 0.508, R:0.0165, max dev:0.0161 


АЛ, 


Figure 2.16 Fit of. æwz(PTFY/ÆTF,, 
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1.50 











The resulting force and moment model in the launchway are expressed as 


CFY = 7. (асвуртву аль үү )А/ 1 Д(СЕҮ, енен) 
CMZ = fal ame (PIE) PIE, ¡JA! L+ (СМ) 


where f, represents a linear fit 
f, represents a third order polynomial fit 


CFY(PTF) for A/L 1.0 





Figure 2.17 Sample of Lateral Force Model Fit A/L = 1.0 
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CMZ(PTF) for A/L 1.0 





Figure 2.18 Sample of Lateral Moment Model Fit A/L = 1.0 


2.2.2 Host Proximity 


With the vehicle near the boundary but not in the launchway, insufficient data 
existed in [11] to predict the effects of relative pitch and clearance from the boundary. 
To develop this model, an inviscid potential flow panel code was resorted to, in order to 
determine the near host interaction effects. Difficulties in the numerical hydrodynamics 
were encountered in modeling due to the blunt forward end ofthe UUV geometry. Initial 
paneled models attempted to reflect the bow accurately, however the large change in 
sectional area vs. length resulted in a large negative pressure at the bow leading to 
excessive axial suction force instead ofthe zero drag condition expected. Modification of 
the model to blend the change in bow cross section over increased length provided a 
model which corresponded well with strip theory and experimental data to a similar 


vehicle [10]. The resulting panel geometry is illustrated in Figure 2.19. 


28 











Figure 2.19 UUV Panel Geometry 


The proximity model is obtained in two parts. First, investigating the effect of relative 
proximity to the boundary both with and without relative pitch including far field 
determination of the Munk moment using potential flow. And second, investigating the 
effect of PTF on force and moment from experimental data. Figure 2.20 shows the 
suction force and moment as a function of gap/diameter, where gap refers to the distance 
from the vehicle molded shell to the boundary. The effect of pitch and gap/diameter on 
the force and moment 1s illustrated for selected cases in Figures 2.21 through 2.23. The 
non linear least square fit is superimposed in Figures 2.21 and Figure 2.22 showing good 
agreement over the range of combinations of pitch and gap/diameter. Figure 2.23 
represents the potential flow Munk moment determined with the vehicle in an unbounded 


fluid domain. 
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CFY and CMZ At Zero Pitch 
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ean: y-a* exp (b* x". 6), R:0.000122, max dev:3.63E-4 736948, b --107.1891 
— — — - eqn: y-a* exp(b*x ), R:0.000383, max dev:&.23E- 4 27-0.484992, Ь--3.00457, с-0.570998 
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gap/dia 


Figure 2.20 Boundary Interaction At Zero Pitch 


CFY vs gap/dia for Pitch From 0° to 20° 


Least Squares Fit (6* erder poly) 
Potential Flow 


gap/dia 


Figure 2.21 Lateral Force Near a Boundary 
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CMZ vs gap/diameter for Pitch of 1? 
0.050 


о Least Squares Fit 
o Potential Flow 
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Figure 2.22 Lateral Moment Near a Boundary 
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Figure 2.23 Open Water Munk Moment 
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The effect of PTF on the force and moment is captured by data obtained in [11]. Figures 
2.24 and 2.25 represent the lateral force and moment effects due to a change in propulsive 
thrust factor. Significant scatter of experimental data was observed and several outlier 
points have been omitted in both curve fits, using the method described previously, to 


arrive at a model which makes physical sense and tends to the proper limit. 


CFY Near Wall 


eqn: y=a*x, R:0.127, max 4еу: 0.173 а=-0.0672913 
d/g = 1.95 

d/g = 31.2 

d/g = 7.8 





Figure 2.24 PTF Effect on Lateral Force 
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CMZ(PT F) 


CMZ(PTF) Near Wall 


eqn: y=a*exp(b* Ux), R:0.00178, max dev:0.00245 
a=4.13743, b=-156.197 





dia/gap 


Figure 2.25 PTF Effect on Lateral Moment 
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Chapter 3 SIMULATION MODEL 


3.1 Coordinate Systems 


The forces acting on an underwater vehicle are normally referred to a non-inertial 
coordinate system fixed in the body, however, motion ofthe vehicle is most easily 
described in an inertial reference frame. For the purposes of this analysis, the inertial 
reference frame is translating at the constant velocity ofthe host, assumed to be along the 
positive X axis, and therefore represents a host fixed reference. The Z axis is directed 
positive downwards and the positive Y axis is to starboard in this frame of reference. To 
obtain the forces and moments in the non-imertial reference it is necessary to transform the 
vehicle motion to this frame. The vector transformation between the reference frames is 
uniquely determined by three angles: > rotation about the X axis (roll), 0 rotation about 
the Y axis (pitch), v rotation about the Z axis (yaw), which describe the angular 
orientation between the two frames and are termed the Euler angles. The order in which 
the transformation about these angular relations is observed is significant, and if taken in 
different sequence can produce a different vector result. In accordance with standard 
convention the transformation is accomplished by ordering the rotations as first yaw ( vy), 
then pitch (0), and finally roll($). The transformation matrix describing the vehicle motion 


in translation in the non-inertial frame, given the inertial frame motion, is expressed by: 


cos( 0) cos( y) cos( 0) sin( v) —sin( Q) 
T, (9, 8, vy) =| —sin(y) cos( g) + sin(g)sin(@)cos(y)  cos(ó)cos(y)-* sin(ó)sin(O)sin(y) sin(ø)cos(0) 
sin( Ø) sin( vw) +cos(ø)cos(w)sin(0) —sin(¢)cos( y)+cos(¢)sin(@)sin(y) cos( ¢)cos( 8) 


This transformation accounts for the rigid body rotation from host fixed reference to body 
fixed reference. The transformation of angular rate quantities follows the same Euler angle 


transformation convention and order of rotations. 
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The angular rate transformation is then: 


1 0 — $11( 0) 
T,(8,0,y)=|0 сов(0) sin(ø)cos(ø) 
0 -sin(d) cos(O)cos(9) 


The total transformation matrix may be represented as: 


T, (é, 0, v) 0 | 
0 7) (9,0, v) 


T(9,0, yw) = | 
With R, representing a vector in the inertial frame consisting of scalar quantities of u 
(surge), v (sway), and w (heave) velocities , and p (roll rate), q (pitch rate) and r (yaw 
rate), or their respective time derivatives; the previously defined Euler angles $, 9, and y 
are used to determine the vector ®, representing the same vector 91,, but in the body 
reference frame changed only in orientation. The relation 

$t - T($,0, v)8 

is used to determine the body relative velocity and acceleration in the body reference 
frame translating with the velocity of the host in the host axial direction. Finally, the body 


velocities are corrected for this constant translation of the body fixed reference frame to 


arrive at the total body relative velocities independent of the host. 
3.2 Equations of Motion and Implementation 


The equations of motion of a rigid body may be expressed for arbitrary origin in 


matrix form as follows: 


шэн 0 0 0 mzg -тУс 7 X 
0 m— Y; 0 Бин 0 mxg —Y, | v 

0 0 Maza тУс —mx ç = 2; 0 w Z 

0 —mz ç — K; тУс iie = E р K 

т? с 0 -тхс – М, и Le Vee EE M 
|-myg тхо №, 0 -1,-М, -1, 1,-М;|4| |м| 
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This expresses the Revised Standard Equations of Motion [5] with the acceleration terms 
isolated to the left hand side. The right hand side consists of the sum of the gyroscopic 
and centripetal acceleration effects as well as forces arising from the fluid hydrodynamic 
effects on the vehicle. These last include forces and moments accounted for by 
maneuvering coefficients including control forces, hydrostatics, crossflow and vortex 
effects. The hydrodynamic forces are calculated in accordance with the standard 


equations of motion. These are described for the hydrodynamic terms as: 


X=X,,q° +X_r? риа +X v? + سر‎ + Хи» + Х ид + Х и“ 


= ۲, |p| + Ра Ра + Бир Кира mur Y uve] t 1214 Байр + | + Frie] + uu + kin 


Му? + и? 





+» 


Z = Z¿uq + Z,uw + Z| + Дур + дьии + ЖЕСІ - Zum w| w| +Z pt Z M 








K=K,up+ K,ur+K,,qr+K,,wp+K,,wr +K,uv+K ; p|p|+ K.u? 
M = M,uq +M „„alal+ M „„wlw|+ M,,vp + Myr + Mau? + M,,uw+M,,, wv? +? 


12712 


М = Мара + Мрир + Мару + Мен“ + Мит + М w+N ur 








ЕБ Мир ar МУ“ + Ny vy ен N „rlr 


The body terms, including coriolis, centripetal, and hydrostatic effects are calculated 


separately. Neglecting cross products of inertia these terms may be expressed in the form: 


Хуу = стнд -vr = xg (q? +r) +ygpq + 2g (pF))-( = B)sin(6) 

Frody = =m(ur — wp + xg pq + 2gqr)+(W — B)cos(6)sin( $) 

©, - -mlvp — ug +xGpr-25(p" +9?) + (1 - B)cos(6)cos( 6) 

Kbody = U yy — 172 Jar — ml ya (vp — q ) +26 (ur — wp )) + (y W — ygB)cos( 9) cos( ó) — (zGW — zgB)cos(6) sin( 4) 
M body UL )pr - m(xg (vp —uq) - zg (wq —-vr)) - (xgW — xgB)cos(6) cos( 9) —(2¿W —z5B)sin(6) 

Мау = (1, — 1 yy) pq ^ m(xg(ur — wp) - ya(wq —vr)) * (xoW —x5B)cos(6)sin(4) + (Yo — ygB)sin(6) 
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The hydrodynamic maneuvering coefficients may be obtained by experiment, or predicted 
by numerous methods [1], [9], [12]. These relations do not include the propulsive thrust 
model forces or control forces in order to prevent redundancy of these terms with the 
interaction force model. 

The right hand side forces on the vehicle may be expressed as a partition of the 
total right hand side in the following fashion: 

FE. PIE ES 
where F,,,,, may be equated to the acceleration and mass product terms , F,, represents the 
forces and moments acting on the body due to open water hydrodynamic, centripetal, 
gyroscopic, and hydrostatic effects. The term F, accounts for the interaction effects due to 
the operation in a restricted fluid e.g. in the launchway and host vicinity. F, accounts for 
the control force and moment required to effect the selected trajectory. F, may be solved 
for directly: 

Fo E, uper Ep 
This relation is adapted for implementation in a six degree of freedom motion simulation 
program for which a general description of the method employed follows. 

The desired trajectory of the vehicle must be known a - priori and specified as 
input to the algorithm. It is important in the launch and recovery of the vehicle to know 
the relative position , velocity, and acceleration of the vehicle with respect to the moving 
host, therefore the trajectory is specified in a host fixed inertial frame. The vehicle 
accelerations and velocities are determined by a coordinate transformation from the 
moving inertial frame to the body fixed frame with the effect of forward speed of the host 
accounted for in the transformation. The total force at each time step 1s determined from 
the vehicle accelerations by solving the product of the mass matrix and body relative 
acceleration vector directly. A word of explanation is warranted with respect to the use of 
a motion simulation program. For the aspect of this thesis, with the specification of a 


trajectory known explicitly over time along with acceleration, and velocity, there is no 
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need to invert the mass matrix and integrate over each time step. The forces and moments 
may be computed directly at each time step allowing the algorithm to run very fast without 
growing errors due to a integration scheme. The hydrodynamic forces and moments 
including the body terms are calculated at each time step without consideration to the host 
presence. A separate subroutine implements the interaction force model previously 
developed, where forces and moments when the vehicle is located within the launchway 
are determined completely from the launchway model and include hydrodynamic and 
interaction effects. Once the vehicle has transitioned to a position outside of the 
launchway the interaction effects alone are determined. This condition reflects only a 
portion of the forces acting upon the vehicle, therefore to arrive at the total right hand side 
force vector when the vehicle is outside of the launchway, the hydrodynamic derivative 
based forces are added to the interaction effects. 

In order to achieve a smooth transition between the experimental data representing 
the total force on the vehicle to the force calculated by the hydrodynamic derivative based 
force vector, it is necessary to employ a smoothing procedure to blend the differences 
between the two. In order to accomplish this, a zone structure is defined with the 


following description: 


Zone 1 Inside the host launchway F(x,y,z) - Hydrodynamic, + Interaction 
Zone2 Within the transition region F(x,y,z)= Hydrodynamic; + Interaction; 
Zone3 Outside of the transition region F(x,y,z) = Hydrodynamicg + Interactione 
Where F(x,y,z) is the force calculated by the force model and subscripts 
e, i, and q refer to experiment, interpolated, and coefficient based 


semi-empirical data respectively. 


Zone 2 is where the smoothing routine operates. This zone is defined by a hemisphere 


with the center located at the launchway exit defined by x¿,yc,Z¿ With the diameter of this 
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4) Å geometry limitation file containing the position information with regard 
to launchway geometry, host geometry and host speed. 
5) A trajectory file which includes vehicle acceleration and velocity scaling 


parameters and trajectory modification commands versus time. 


A linear interpolation scheme is used to determine the body position, velocity and 
acceleration changes over time from the commanded trajectory inputs. A partial list of the 
subroutines implementing the trajectory, force model and coordinate transformation is 


included as Appendix A. 


40 








Chapter 4 SIMULATION RESULTS 


4.1 Trajectory Runs 


Four different trajectories were selected for evaluation of the simulation. A copy 
of the input files is included in Appendix B. The trajectories were selected to test the 
blending between the force and moment models between Zone 1 and Zone 3. The 
trajectories selected were not intended to reflect in all cases the trajectory that would be 


optimal or would represent the real vehicle desired course but would illustrate the model. 
4.2 Straight Launch 


The first trajectory simulated was a straight ahead motion of the vehicle. Using an 


exponential time function acceleration described by: 
x(t) = 2(1- e" ) 


2 
Xpoa, (1) = (4*0. en 


) 
The vehicle is directed out ofthe launchway and away from the host at a constant heading 
angle. The trajectory in host reference is illustrated in Figure 4.1. Figure 4.2 through 4.7 


show the resulting accelerations, velocities, forces and moments for this maneuver in the 


vehicle reference frame. 
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Figure 4.1 Vehicle Trajectory 


Traj 1 Body Relative Velocities 
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Figure 4.2 Velocity Profile 
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Traj 1 Vehicle Relative Acceleration 
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Figure 4.3 Acceleration Profile 
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Figure 4.4 Hydrodynamic Forces in Open Water 


43 





Force [Ib] 


Force [lb] 


Traj 1 Interaction and Hydrodynamic Forces 
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Figure 4.5 Hydrodynamic Forces Including Vehicle Interactions 
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Figure 4.6 Total Force 
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Traj 1 Control Force 
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Figure 4.7 Control Forces 


Although unrealistic as a desired maneuver due to the large sustained crossflow drag 
resulting from the constrained yaw angle, this trajectory is illustrative of the interaction 
effects and transition from a region where interaction effects play a role to the region in 
the far field where the interaction effects are small. As such it is useful as a baseline 
condition. The location of the vehicle relative to the zones previously defined is indicated 
in the figures, where pertinent, by a zone indicator line which takes a value of 1, 2, or 3 
respective to the appropriate zone. 

The interaction effects of the launchway and host vehicle are apparent in the figures. As 
the vehicle nose passes the launchway exit (occurring very early in the trajectory at — 4 
sec) the lateral force can be seen to decrease as the nose proceeds into the transition zone. 
The effect on the lateral moment indicates an almost impulsive effect turning the moment 
from a negative value when the vehicle is fully within the launchway to a large positive 


moment as the vehicle fully exits. The magnitude of this moment depends upon the 
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interaction due to the boundary as well as the effect of PTF on the accelerated flow 
between the vehicle and the interacting boundary. As the boundary recedes, the moment 
decreases to the open water Munk moment due to the imposed angle of attack of the 


vehicle in this trajectory simulation. 


4.3 Launch with Yaw Maneuver 


The second trajectory for simulation yaws the vehicle to a zero relative heading 
angle after the vehicle has cleared the launchway. This represents a more realistic 
maneuver intended to investigate the effect of maneuvering to remove crossflow. The 
vehicle is yawed to a zero heading angle at time (t = 25 secj and proceeds along a path 
parallel to the host for the rest of the simulation. This trajectory is illustrated in Figure 


4.8. Figures 4.9 through 4.14 illustrate the velocities, forces and moments. 
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Figure 4.8 Vehicle Trajectory 
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Traj 2 Body Relative Velocities 
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Figure 4.9 Velocity Profile 
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Figure 4.10 Acceleration Profile 
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Figure 4.11 Hydrodynamic Forces in Open Water 
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Figure 4.12 Hydrodynamic Forces Including Vehicle Interactions 
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Figure 4.13 Total Force 
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Figure 4.14 Control Force 
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The interaction effects near the host are apparent in the lateral force and moment as they 
were in the straight launch trajectory. In addition, acceleration effects during the 
maneuver can be seen to affect the required control forces, specifically in this instance the 


contribution of V; to the lateral moment. 


4.4 Launch with Yaw and Deceleration Maneuver 


The third trajectory selected for simulation illustrates a maneuver intended to 
reflect a realistic vehicle maneuver. In this simulation the vehicle transitions out ofthe 
launchway and subsequently decelerates to the desired speed once clear of the host. 
Figure 4.15 shows the trajectory history. Figures 4.16 through 4.21 illustrate the 


velocities, forces and moments. 
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Figure 4.15 Vehicle Trajectory 
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Figure 4.17 Acceleration Profile 
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Figure 4.18 Hydrodynamic Forces in Open Water 
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Figure 4.19 Hydrodynamic Forces Including Vehicle Interactions 
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Figure 4.20 Total Force 


Traj 3 Control Force 





[ат 3) лошој 


Ч 


-80 


-160 


100 





80 


20 


Time [sec] 


Figure 4.21 Control Force 
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Again as in the previous trajectory the significance of the acceleration effects are evident 
in the axial and lateral body forces during deceleration. The lateral moment peak occurs 
just after the vehicle has cleared the launchway and experiences the maximum interaction 


effect of the planar boundary in proximity to the vehicle in this configuration. 


4.5 Straight Launch At Zero Host Speed 


The fourth trajectory selected for simulation illustrates a launch from a stationary 
host. In this simulation the vehicle transitions out ofthe launchway maintaining a straight 
course and speed without onset flow. Figure 4.22 shows the trajectory history. Figures 


4.23 through 4.28 illustrate the velocities, forces and moments. 
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Figure 4.22 Vehicle Trajectory 
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Figure 4.23 Velocity Profile 
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Figure 4.24 Acceleration Profile 
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Figure 4.25 Hydrodynamic Forces in Open Water 
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Figure 4.26 Hydrodynamic Forces Including Vehicle Interactions 
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Figure 4.27 Total Force 
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Figure 4.28 Control Force 
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It is of interest to note the differences between this launch and the previous launches, with 
forward host speed. In this launch the interaction effects are clearly seen; the lateral force 
changing sign from an outward force when the vehicle is mostly within the launchway to a 
inward suction force once clear. The lateral moment also changes sign from a bow in 
moment when the vehicle is mostly within the launchway to an outward moment due to 
boundary interaction and suction region acting on the aft portion of the vehicle. The axial 
force reflects the acceleration to constant velocity with corresponding constant axial 
hydrodynamic drag, while the lateral moment approaches asymptotically to a value 
dependent upon PTF. The effect of forward host speed has a marked effect on the force 
and moment on the vehicle, from commencement of launch with the vehicle fully within 
the launchway, to its transition to open water. With a 4 feet per second onset flow 
velocity due to the host forward speed, the lateral force does not achieve sign reversal, 
and an outward force dominates throughout the maneuver sequence. The effect on 
moment is pronounced as well, reflecting a twenty fold increase in peak moment over the 
zero forward speed case. 

The small perturbation of the lateral moment observed as the vehicle becomes fully 
emerged from the launchway is qualitatively similar to the moment profile results 
described in [2] for a high speed torpedo ejection from a submarine launchway. This 
comparison is not intended to validate the results here but is offered for completeness. 
Due to the disimilar nature of the ejection dynamics compared to vehicle swimout, a 


quantitative comparison was not conducted. 
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Chapter 5 CONCLUSION 


5.1 Summary 


From the comparison of trajectory runs it is clear that the force and moment model 
implemented into a six degree of freedom motion algorithm can solve for control forces 
versus time for any desired trajectory. The mathematical model successfully makes non 
dimensional and integrates experimental data with potential flow panel code results, and 
with a parametric based non linear mathematical model, to span the trajectory regime from 
launchway to open water. Host vessel speed effects are incorporated. The smoothing 
technique accounts for the inevitable mismatch between the experimentally derived 
interaction model to the parametric, coefficient based, hydrodynamic forces in open water 
and allows for geometrically similar models to be tested using the existing experimental 


results. 


Based upon trajectories with changing accelerations, acceleration effects add 
significant information to the force simulation. Vehicle motion gain evaluation is enabled 


through the trajectory input file. 


The magnitudes of the required control forces are broken down in such a way so 
as to enable the contributing interaction forces, open water hydrodynamic forces, or body 


mass forces, to be compared and their dependence upon trajectory evaluated. 


5.2 Recommendations for Future Work 


The control force output from the model represents, in a sense, the magnitude of 


control actuated force to maintain the selected trajectory. A shortcoming in the model 
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exists in that it ignores interaction effects between the control surface or thruster, the 
vehicle, and proximity boundary. The fidelity of the simulation would be improved by 


incorporating control effects on the uncontrolled vehicle interaction force model. 


For a vehicle with a different propulsor, a means of effectively adapting the thrust 


model and resulting thrust, moment and torque effects would be useful. 


The basis for a maneuvering workbench exists which can be used to assist in 
control appendage sizing, control system design, and free swimming motion simulation in 


one simulation environment. 


The accuracy of the simulation and its representation of the real world can only be 
determined by comparing its output to full scale tests. These would indicate those areas 
captured by the model as well as areas requiring improvement. Additional scale testing 
should be conducted to provide a more robust mathematical model, and should 


incorporate thruster, propeller, and launchway off - axis effects. 
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Appendix A. Source Code Listing for Trajectory Simulation and Force Model 
This appendix includes source code for the implementation of the trajectory, and force 


model. Also included is a portion of the Equation of Motion (EOM) program used to 
calculate the hydrodynamic open water forces. 
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С Ad ode ode ode old ad ole ole ole ode ode ade al ale ale ale ale ale ode ode ode ade ale ale ale ale ode ode ade ale ale ale ale ade ade ale ale ale ade ade ae ade ale ode ade ade ade ale ale ale k k kk k k ad ad k k al k kk k 


SUBROUTINE TRJ 


IMPLICIT NONE 

REAL*8 CTRN(6,6),RHS(6),EDOT(6),E(6),PDOT(6) 

REAL*8 P(6),DT, TA,PDOTO,PDDOTO 

REAL*8 DS4,DC4,DS5,DC5,DS6,DC6,DT5,DS5DC6,DS5DS6 
REAL*8 PDDOT(6),RHSTOT(6),RHSC(6) 

REAL*8 A(6,6), WORK(6,6),ACOF(6),BCOF(6) 

REAL*8 ACP(3),BCP(3),RHSFM(6),RTD,DTR 

REAL*8 XC,YC,ZC,XMIN,XMAX, YMIN, YMAX,ZMIN,ZMAX,D2,HRAD 
REAL*8 VRAD, UHOST,DIA,RHOC,ATLMAX,VLEN 
REAL*8 FX,FY,FZ,MX,MY,MZ,AEDOTI,RHS1,PTF 
REAL*8 ATL, VEX(2), TEMP1,TEMP2 

REAL*8 P1A,P1T,P2A,P2T,P3A,P3T 

REAL*8 P4A,P4T,P5A,P5T,P6A,P6T 

REAL*8 PIDOT,P2DOT,P3DOT 

REAL*8 P4DOT,P5DOT,P6DOT 

REAL*8 PIDDOT,P2DDOT,P3DDOT 

REAL*8 PADDOT,P5DDOT,P6DDOT 

INTEGER*2 LJ,J,K,JJJ,ILKODE,JR,JJR,JP,JJP,JYA,JJYA 
INTEGER®*2 JX,JIX,JY,JJY,JZ,JIZ 

CHARACTER*15 LABEL 


REAL*8 DD(1800) 

COMMON /COMDD/ DD 

COMMON/FGEOM/ XC, YC,ZC,XMIN, XMAX,YMIN, YMAX,ZMIN,ZMAX, 
1 D2, HRAD, VRAD, UHOST,DIA,RHOC,ATLMAX, VLEN 
1 ACP,ACOF,BCP,BCOF 

SAVE JJIJRJIRJPJIP,JYAJJYA,JX,JIX,JY,JIY,JZ,JIZ 

DATA J,J],JR,JJR,JP,JJP,JY A,JJY A /9,9,2,0,2,0,2,0/ 

DATA JX,JJX JY JJY JZ_JIZ /2,0,2,0,2,0/ 


EQUIVALENCE (DD( 1), DT) 

EQUIVALENCE (DD( 11), JJJ),DD( 4, RTD) 
EQUIVALENCE (DD( 10), — TA,(DD( 5, DTR) 
EQUIVALENCE (DD( 31),EDOT(1)) 
EQUIVALENCE (DD( 55), E(1)) 
EQUIVALENCE (DD( 61), Р(1)) 
EQUIVALENCE (DD( 71), RHS(1)) 
EQUIVALENCE (DD( 77),PDOT(1)) 
EQUIVALENCE (DD(107),A(1,1)) 

EQUIVALENCE (DD(83),PDDOT(1)) 
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DIMENSION P1A(10),P1T(10),P2A(10),P2T(10),P3A(10),P3T(10) 
DIMENSION P4A(10),P4T(10),P5A(10),P5T(10),P6A(10),P6T(10) 


IF (TA .EQ. 0.) THEN 
C READ IN THE GEOMETRY CONSTRAINTS 


WRITE(6,*) ' OPENING AND READING GEOMETRY LIMITATION FILE’ 
OPEN(11,FILE='GEOM.LIM',STATUS='OLD’) 

READ(11,*) XC,YC,ZC 

READ(11,*) XMIN, YMIN,ZMIN 

READ(11,*) XMAX, YMAX,ZMAX 

READ(11,*) D2, HRAD, VRAD, VLEN 

READ(11,*) UHOST,DIA,RHOC,ATLMAX 

CLOSE (11) 


C READ IN FORCE MODEL COEFFICIENTS 


WRITE(6,*) ' OPENING AND READING LEAST SQUARES COEFFICIENT 
FILE 

OPEN(12,FILE=Z2COEF.EXP', STATUS=OLD') 

READ (12, *) (АСР(П), П = 1,3) 

ВЕАР” (12,*) (ВСРП), П = 1,3) 

READ (12,*) (АСОЕ(П), П = 1,6) 

READ (12,*) (ВСОЕ(П), П = 1,6) 

CLOSE (12) 


C READ IN TRAJECTORY MODIFICATIONS 


WRITE(6,*) ' OPENING AND READING TRAJECTORY COMMAND FILE' 
OPEN(14,FILE='TRAJ.CM', STATUS='OLD’) 
READ (14,*) (УЕХ(П), П = 1,2) 
READ (14,*) 
READ (14,*) 
DO 20 I = 1,10 
КЕАГ(14,%) РІА(1),Р1Т(1),Р2А(1),Р2Т(1),РЗА(1),РЗТ(1) 
20 CONTINUE 
READ(14,*) 
READ(14,*) 
DO 30 I = 1,10 
READ(14,*) РААП),РАТ(1) РБА(1) Р5Т(1) РбБА(1),РбТ(1) 
P4A(I) = P4A(I) * DTR 
PSA(I) = PSA(I) * DTR 
P6A(I) = РбА(1) * DTR 
30 CONTINUE 
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CLOSE (14) 


C COMPLETE MODEL ENTRY ` 
ELSE 
ENDIF 


C USE STANDARD INITIAL VELOCITY PROFILE 
C MODIFY TRAJECTORY BY SCALING OF U,V, AND W RELATIVE 
C TO THE HOST FOLLOWED BY ROLL,PITCH,YAW RATE COMMANDS 


IF (TA .LT. PIT(JX)) THEN 

ID JX-1 

PIDDOT - (PIAQX) - PIAQJX)/(P1T(UX) - P1T(JX)) 
PIDOT = PIDOT + PIDDOT*DT 

ELSE 

PIDDOT = (PIA(JX) - PIA(JJX))/(P1T(JX) - PIT(JJX)) 
PIDOT = PIDOT + PIDDOT*DT 

ЇХ-1Х-1 

ENDIF 


C  WRITE(6,*) TA,PIDOT,PIDDOT 


IF (TA .LT. P2T(JY)) THEN 

WY=JY-1 

P2DDOT = (P2A(JY) - P2A(JTY))/(P2T(SY) - P2T(ITY)) 
P2DOT = P2DOT + P2DDOT*DT 

ELSE 

P2DDOT = (P2A(JY) - P2A(JY))(P2TQY) - P2T(JY)) 
P2DOT = P2DOT + P2DDOT*DT 

ТҮ-1Ү-1 

ENDIF 


IF (TA .LT. P3T(JZ)) THEN 

JZ =9Z- 1 

P3DDOT = (P3A(JZ) - P3A(JIZ))(P3T(JZ) - P3T(JJZ)) 
P3DOT = P3DOT + P3DDOT*DT 

ELSE 

P3DDOT = (P3A(JZ) - P3A(JJZ))/(P3T(JZ) - P3T(JJZ)) 
P3DOT = P3DOT + P3DDOT*DT 

JZ=JZ +1 

ENDIF 


IF (TA .LT. P4T(JR)) THEN 


JR —SR- 1 
P4DDOT = (P4A(JR) - PAA(JJR)/(PAT(JR) - PAT(JJR)) 
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P4DOT = P4DOT + PADDOT*DT 

ELSE 

P4DDOT = (P4A(JR) - P4A(JJR))/(P4TUR) - PAT(JIR)) 
P4DOT = P4DOT + PADDOT*DT 

JR = ЈЕ + 1 

ENDIF 


IF (TA .LT. PST(JP)) THEN 

JJP - JP - 1 

P5DDOT = (P5A(JP) - PSA(JJP))(PST(JP) - PST(JJP)) 
P5DOT - P5DOT -- PSDDOT*DT 

ELSE 

PSDDOT = (P5A(JP) - PSA(JJP))/(PST(JP) - PST(JJP)) 
PSDOT = P5DOT + PSDDOT*DT 

JP=JP+1 

ENDIF 


IF (TA .LT. P6T(JYA)) THEN 

ЈУА = ЈУА - 1 

P6DDOT - (P6A(JYA) - PeEA(JJYA))(P6T(JYA) - PGT(JJYA)) 
P6DOT = P6DOT + P6DDOT*DT 

ELSE 

P6DDOT = (P6A(JYA) - PEA(JYA))(P6T(JYA) - PGT(JJYA)) 
P6DOT = P6DOT + P6DDOT*DT 

JYA =JYA +1 

ENDIF 


C ESTABLISH THE BODY POSITION, VELOCITY, AND ACCELERATION 
C OVER TIME IN INERTIAL REFERENCE FRAME COORDINATES 
C RELATIVE TO THE HOST 


PDOTO = VEX(1)*(1.- DEXP(-(VEX(2))*(TA**2))) 
PDOT(1) = PDOTO*DCOS(P(6)) + PIDOT 
PDOT(2) = PDOTO*DSIN(P(6)) + P2DOT 
PDOT(3) = PDOTO*DTAN(P(5)) + P3DOT 
PDOT(4) = PADOT 
PDOT(5) - PSDOT 
PDOT(6) = P6DOT 


PDDOTO - VEX(1)*(2.* VEXQ)*TA*DEXP(-(VEX(2))*(TA**2))) 
PDDOT(1) = PDDOTO*DCOS(P(6)) + PIDDOT 

PDDOT(2) = PDDOTO*DSIN(P(6)) + P2DDOT 

PDDOT(3) = PDDOTO*DTAN(P(5)) + P3DDOT 

PDDOT(4) = PADDOT 

PDDOT(5) = PSDDOT 
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PDDOT(6) = PODDOT 


P(1) = P(1) + PDOT(1)*DT + PDDOT(1)*DT**2 
P(2) = P(2) + PDOT(2)*DT + PDDOT(2)*DT**2 
P(3) = P(3) + PDOT(3)*DT + PDDOT(3)*DT**2 
P(4) = P(4) + PDOT(4)*DT + PDDOT(4)*DT**2 
P(5) = P(5) + PDOT(5)*DT + PDDOT(5)*DT**2 
P(6) = P(6) + PDOT(6)*DT + PDDOT(6)*DT**2 


C CALCULATE COORDINATE TRANSFORM FROM INERTIA FRAME (HOST) 
C TO VEHICLE FRAME 

C TRANSPOSE OF THE SNAME 1950 CONVENTION (IAW ABKOWICZ) 

C ORDER IS YAW, PITCH, ROLL 


DS4 = DSIN(P(4)) 
DC4 = DCOS(P(4)) 
085 = DSIN(P(5)) 
DC5 = DCOS(P(5)) 
DS6 = DSIN(P(6)) 
DC6 = DCOS(P(6)) 


CTRN(1,1)= DC5*DC6 

CTRN(1,2)= DC5*DS6 
CTRN(1,3)=-DS5 

CTRN(2,1) = -DS6*DC4 + DS4*DS5*DC6 
CTRN(2,2) = DC4*DC6 + DS4*DS5*DS6 
CTRN(2,3)= DS4*DC5 

CTRN(3,1)= DS4*DS6 +DC4*DC6*DS5 
CTRN(3,2) = -DS4*DC6 + DC4*DS5*DS6 
CTRN(3,3) = DC4*DC5 


CTRN(4,4) = 1.0 
CTRN(4,5) = 0.0 
CTRN(4,6) = -DS5 
CTRN(5,4) = 0.0 
CTRN(5,5) = DC4 
CTRN(5,6) = DS4*DC5 
CTRN(6,4) = 0.0 
CTRN(6,5) = -DS4 
CTRN(6,6)= DC5*DC4 


С XFORM INERTIA RELATIVE VELOCITIES AND ACCELERATIONS 
C TO VEHICLE RELATIVE 

С CTRN (6X6) * PDOT (6X1) = E (6X1) 

C CTRN (6X6) * PDDOT (6X1) = ЕРОТ (6X1) 
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C INITIALIZE THE E AND EDOT MATRICE TO ZERO AT EACH TIME STEP 


DO 110 I=1,6 

E(D = 0.0 

EDOT(D = 0.0 
110 CONTINUE 


DO 120 I=1,6 

DO 130 K=1,6 

E(D = E(I) + PDOT(K)*CTRN(LK) 

EDOT(I) = EDOT(I) + PDDOT(K)*CTRN(LK) 
130 CONTINUE 
120 CONTINUE 


C CONVERT BODY VELOCITIES TO ADD IN HOST VESSEL SPEED EFFECTS 
C HOST SPEED IN (-) P(1) DIRECTION ONLY. ADD AS A POSITIVE SCALAR 


E(1) = E(1) + UHOST*DCOS(P(6))*DCOS(P(5)) 
E(2) - EQ) - UHOST*DSIN(P(6)) 
E(3) - E(3) + UHOST*DCOS(P(6))*DSIN(P(5)) 


C SOLVE FOR THE RIGHT HAND SIDE VECTOR WITH TOTAL FORCE 
C EDOT (6X1) * LHS (6X6) = RHS (6X1) 


C INITIALIZE THE RHSTOT TO ZERO AT EACH TIME STEP 
DO 140 I=1,6 
RHSTOT(I) = 0.0 
140 CONTINUE 


IF (JJ .EQ. 9) THEN 
WRITE(25,178) (acp(I),I=1,3) 


LABEL = 'CTRN' 
CALL MPRINTI(CTRN,6,6,9, LABEL) 


LABEL - 'LHS' 

CALL MPRINT1(A,6,6,9, LABEL) 
ENDIF 

JJ=JJ+1 


DO 200 I=1,6 


DO 220 K=1,6 
RHSTOT(I) - RHSTOT(I) -A(LK)*EDOT(K) 
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220 CONTINUE 
200 CONTINUE 
AEDOTI - RHSTOT(1) 


C CALCULATE THE RIGHT HAND SIDE OF THE HYDRODYNAMIC 
EQUATIONS 
C WITH NO CONTROL FORCE OR INTERACTION FORCE 


C SET CROSSFLOW SWITCH JJJ TO 1 


JJJ=1 
CALL HYDRO 


C CALCULATE THE INTERACTION FORCES IN FORCE MODEL ROUTINE 
CALL FM (FX,FY,FZ,MX,MY,MZ,AEDOTI,PTF,ATL,KODE) 


C CALCULATE THE CONTROL FORCES TO EFFECT THE 
C CHOSEN TRAJECTORY. ADD IN HYDRO FORCE 


RHSFM(1) = FX 
RHSFM(2) = FY 
RHSFM(3) = FZ 
RHSFM(4) = MX 
ЕНЗЕМ(5) = МУ 
RHSFM(6) = MZ 


DO 240 1-1,6 
RHSC(I) = RHSTOT(I) - RHSFM(I) 
240 CONTINUE 


C PRINT OUTPUT EVERY 10 DELTA T 


J= JH 
IF (J .EQ. 10.) THEN 
J=0 


WRITE(21,176) TA, (RHSTOT(I),I=1,6) 

WRITE(22,178) TA, KODE, (RHS(I),I=1,6) 

WRITE(23,178) ТА, KODE, (RHSFM(I),I=1,6) 
WRITE(24,177) TA, KODE, (RHSC(I),I=1,6) 

WRITE(25,176) TA,(P(1),I=1,3) (RTD*P(1),I=4,6) 
WRITE(26,176) TA,(PDOT(I),I=1,3),(RTD*PDOT(I),I-4,6) 
WRITE(27,176) TA, (PDDOT(I),I=1,3),(RTD*PDDOT(I),I=4,6) 
WRITE(28, 176) TA,(E(I),I=1,3),(RTD*E(D),I=4,6) 
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WRITE(29,176) TA, (EDOT(I), I 1,3) (RTD*EDOT(I), I=4,6) 
176 FORMAT(ES.1,1X,6(E11.4E3,1X)) 
177 FORMAT(F6.1,1X,12,1X,6(E11.3E3,1X)) 
178 FORMAT(F6. 1,1X,2,1X,6(E9.3,1X)) 

ENDIF 


RETURN 
END 
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C >x k k k < > 2F 2F 2F K K E E E E E K ok a od ol a od ol ol 2k 2k 2 2 K K K de 2K 2K 2F 2K 2K 2K 2K 2 ol R K E R K K R K K K R K K K K E KR K K de k k kK k k 


SUBROUTINE FM (FX,FY,FZ,MX,MY,MZ, 
1 AEDOTI,PTF,ATL,KODE) 
C SUBTOUTINE FORCE MODEL -- VEHICLE SPECIFIC FORCE AND MOMENT 
MODELS ARE PLACED HERE 


CENTER OF X,Y,Z, COORDINATE SYSTEM IS ON HOST AXIS 

IN THIS VERSION, IT IS PRESUMED THAT A NON-CRASH TRAJECTORY 
SPECIFIED 

MODEL IS VALID FOR THE STARBOARD HALF PLANE 

(ALPHA CALCULATION) 


to C © Су © 


REAL*8 MX,MY,MZ,MYC,MZC,MS,MTOT 

REAL*8 P(6),E(6),XC, YC,ZC,XMIN, YMIN,ZMIN 

REAL*8 XMAX,YMAX,ZMAX,D2,HRAD, VRAD,UHOST,DIA 
REAL*8 RHOC,ATLMAX,VLEN,ACP(3),ACOF(6),BCP(3) 
REAL*8 BCOF(6),FX,FY,FZ,DOGS,PTF,PTFS,PTFC 

REAL*8 UZ1LF,UZ1LM,UZ1VF,UZ1VM,FACC,FACS,GAPS 
REAL*8 FYC,FZC,GAP, DOG, PT,SG,FTOT,POSRAD 

REAL*8 UZ2LF,UZ2LM,X,Y,Z,ATL,ATLINV,AXDRAG 
REAL*8 U, V, W,PHI, THETA, PSI, ALPHA, FAREA,C1,C2,C4,C5 
REAL*8 D1,RAT,XS, YS,ZS,AEDOT1,RHS1,TA 

REAL*8 RHS(6) 

REAL*8 DD(1800) 


REAL*4 FOO,FOO1,FOO2,PTF2 


INTEGER KODE 

COMMON /COMDD/ DD 

COMMON /FGEOM/ XC, YC,ZC,XMIN,XMAX, YMIN, YMAX,ZMIN,ZMAX, 
1 D2,HRAD, 

VRAD,UHOST,DIA,RHOC,ATLMAX, VLEN,ACP,ACOF,BCP,BCOF 

EQUIVALENCE (DD( 55), E(1)),(DD( 61), Р(1)) 

EQUIVALENCE (DD( 10), TA)(DD( 71), RHS(1)) 


PI = 3.14159 
DPR = 180.0/pi 
X = P(1) 
Y = P(2) 
Z = PG) 
0-Е(1) 
У = Е(2) 
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W = E(3) 

PHI = P(4) 

ТНЕТА -Р(5) 

PSI = P(6) 

ALPHA = ATAN2(Z Y) 

ЕХ = 0.0 

FY = 0.0 

FZ = 0.0 

MX = 0.0 

MY = 0.0 

MZ = 0.0 

FAREA = PI/4.*DIA**2 

C1 = .5*RHOC*(UHOST**2)*FAREA 
C2 = C1*VLEN 

C4 = .5*(RHOC*U**2)*FAREA 
C5 = C4*VLEN 

C6 = C5/12. 


C IN THIS SUBROUTINE, THE HOST IS TREATED AS A CYLINDER OF RADIUS 
HRAD. 
C POSRAD IS RADIUS FROM HOST AXIS TO SMALL VEHICLE ORIGIN. 


POSRAD = SQRT(Y*Y + Z*Z) 
D1 = SQRT( (X-XC)**2 + (Y-YC)**2 + (Z-ZC)**2 ) 


IF ((Z .LE. ZC) .AND. (Y .LE. YMAX) .AND. (Y .GE. YMIN) .AND. 
I (X.LE. XMAX) .AND. (X .GE. XMIN)) THEN 


ЕЕЕ ЕЕЕ Е ЕЕЕ ЕЕЕ КЕ ОРЕ ЕЕЕ ОЕ ЕЕЕ КОЕ ЕЕЕ ЕЕЕ k k Kk k k k R R K R kR K 


С VEHICLE IS IN THE LAUNCHWAY - ZONE 1 Ч 


NEL EE Eder ЕЕ РЕН ЕЕ ЖЕ ЕЖЕ ЕЕ K Aek ek k kk k K k E k KE EEEk kk Kk REKKER KEKE 


ATL = ATLMAX-(SQRT( (X-XC)**2 + (Y-YC)**2 + (Z-ZC)**2))/VLEN 
ATLINV = ATL - ATLMAX 

AXDRAG = (.637341 +.853415*(ATLINV)**3)*C1 

IF (AXDRAG It. 0.0) AXDRAG = 0.0 

FTOT = AEDOTI + AXDRAG 

PTF = .5*(FTOT)/C1 


FX = -AXDRAG 


UZILF = ((-2.04*ATL**3+5.59*ATL**2-4.49*ATL+1.13)*PTF + 
I (1.40657-.55068*ATL)) 


КУ - UZILF*CI 
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UZILM = -((.153*ATL**3 -.198*ATL**2 - .516*ATL +.508)*PTE) 
MZ = UZILM*C2 


KODE = 1 


ELSEIF ( (Z .LE. ZC) .AND. (Y .GE. YMAX) .AND. (X .GE. XMAX) 
I | .AND. (DI .LT. D2)) THEN 


ЕК ЖЖ ЖК ok ok ok ЕЕ ЖЖ ЖЕ ЖЕ ЖЖ ЖЕ ИК РЕЖЕ ЕЕК ЕС КЕ КЖ ЖУК 


С VEHICLE IS IN TRANSITION ZONE - HEMI-SPHERE ZONE 2  * 
С * 


E ek KNR Е eoe ee ojo ee ek 


IF (D1 .LT. 0.001) D1 = 0.001 
RAT = D2/D1 

XS = XC + RAT*(X-XC) 

YS = YC + RAT*(Y-YC) 

ZS = ZC + RAT*(Z-ZC) 

FACC = (D2-D1)/D2 

FACS = D1/D2 

GAPS = SQRT(YS*YS + ZS*ZS) - HRAD - VRAD 
ЇЕ (GAPS .lt. 0.001) GAPS = 0.001 

DOGS = DIA/GAPS 

ATL = ATLMAX 

ATLINV = 0.0 

AXDRAG = (.637341)*C1 

FXC = AEDOTI + AXDRAG 


C AT LAUNCHER EXIT 


РТЕС = .5*(ЕХС)/С1 


pe MT) 


UZ1LF = ((-2.04* ATL**3+5.59* ATL**2-4.49*ATL+1. 13 )*PTFC + 
1 (1.40657-.55068*ATL)) 


UZILM = -((.153*ATL**3 -.198*ATL**2 - .516*ATL +.508)*PTFC) 


МАС = UZILM*C2 
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C PTF EFFECT AT SPHERE INTERSECTION POINT 
FXS = -RHS(1) + AEDOTI 
PTFS = .5*(FXS)/C4 
UZ2LF = -(0.004*DOGS*PTFS) 
UZ2LM = (4.13743*DEXP(-156.197*(1./DOGS)))*PTFS 
Z2LMPTF = UZ2LM*C5 
C ADD WALL RELATIVE PITCH AND GAP EFFECT 
C APPLY EFFECTIVE PITCH ANGLE DUE TO LAUNCHWAY GEOMETRY 
PTEFF = ((-THETA*SIN(ALPHA) + PSI*COS(ALPHA)) * DPR) 
PT = PTEFF* (XS-XC)/(2.5*VLEN) 
CDBG  WRITE(6,*) PTEFF,PT 
SG = SQRT(DOGS) 
UZ2LF = UZ2LF 
1 -((ACP(1)*DEXP(ACP(2)*(1./DOGS)**ACP(3)))+ 
1 (ACOF(1)*PT*SG+ACOF(2)*PT*SG*SG+ACOF(3)*PT*SG**3 + 
1 ACOF(4)*PT*PT*SG+ACOF(5)*PT*PT*SG*SG+ACOF(6)*PT*PT*SG**3)) 
FS = UZ2LF*C4 
UZ2LM = 
1 (BCOF(1)*PT*SG+BCOF(2)*PT*SG*SG+BCOF(3)*PT*SG**3 + 
1 BCOF(4)*PT*PT*SG+BCOF(5)*PT*PT*SG*SG+BCOF(6)*PT*PT*SG**3) 
MS = UZ2LM*C5+Z2LMPTF 
FX = FACC*(-FXC) + FACS*(-FXS) 


FY = FACC*FYC - FACS*FS*COS(ALPHA) + FACS*RHS(2) 
FZ = FACC*FZC - FACS*FS*SIN(ALPHA) + FACS*RHS(3) 


MY = FACC*MYC - FACS*MS*SIN(ALPHA) + FACS*RHS(5) 
MZ = FACC*MZC + FACS*MS*COS(ALPHA) + FACS*RHS(6) 
PTF = FACC*PTFC + FACS*PTFS 


KODE = 3 
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С 
ELSE 


С K >K K K >K KHK FF KH K > K K K K K K >K K K K >K >K > K K K K > > K K K > K >< K >K > K K K K K K K K K K K K K K 


C ZONE IS ZONE 2 = 


Ce ae ak ala ЖЕЖЕЖЕЖЕЖЕЖЖЖЖЖЖЖЖЖЖЖ 


GAP = SQRT(Y*Y + Z*Z) - HRAD - VRAD 
IF (GAP .LT. 0.001) GAP = 0.001 
DOG = DIA/GAP 
FX = -RHS(1) + AEDOT1 
C PTF EFFECT 
PTF = .5*(-RHS(1) + AEDOT1)/C4 
UZ2LF = -(0.004*DOG*PTF) 
UZ2LM - (4.13743*DEXP(-156.197*(1./DOG)))*PTF 
Z2LMPTF = UZ2LM*C5 
C ADD WALL RELATIVE PITCH AND GAP EFFECT 


PT = (-THETA*SIN(ALPHA) + PSI*COS(ALPHA)) * DPR 
SG = SQRT(DOG) 


UZ2LF = UZ2LF 
1 -((ACP(1)*DEXP(ACP(2)*(1./DOG)**ACP(3)))+ 
1 (ACOF(1)*PT*SG+ACOF(2)*PT* SG*SG+ACOF(3)*PT*SG**3 + 
1 ACOF(4)*PT*PT*SG+ACOF(5)*PT*PT*SG*SG+ACOF(6)*PT*PT*SG**3)) 
FTOT = UZ2LF*C4 
UZ2LM = 
1 (BCOF(1)*PT*SG+BCOF(2)*PT*SG*SG+BCOF(3)*PT*SG**3 + 
1 BCOF(4)*PT*PT*SG+BCOF(5)*PT*PT*SG*SG+BCOF(6)*PT*PT*SG**3) 
MTOT = UZ2LM*C5 + Z2LMPTF 
FX = - FX 
FY = - FTOT*COS(ALPHA) + RHS(2) 


FZ = - FTOT*SIN(ALPHA) + RHS(3) 
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MY = - MTOT*SIN(ALPHA) + RHS(5) 
MZ= MTOT*COS(ALPHA) + RHS(6) 
КОРЕ = 2 

ENDIF 


RETURN 
END 
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С A ode ode ke ade de ole ode ole ode ole ade ole ade ole ade ode ole ade ole ade ole ade ode ole ade ale ade ale ae de ale ade ale ale ade de ole ale ade ale ade ade ole ale de ole ode ole ae ade ole ade ol ae ol ae od ole ol ale de k k k k k k k k 


SUBROUTINE HYDRO 
С жжжжжжжжж жж Жжжж жж жож жж Ж ЖЖ ЖЖ ЖЖЖ ЖЖ ЖЖЖЖ ЖЖЖ ЖЖ ЖЖЖ ЖЖ ЖЖЖ ЖЖ ЖЖ ЖЖ ЖОЖ ЖЖЖ ЖЖ ЖЖ ЖЖЖЖ Ж 
C DTRC 2510 HYDRODYNAMIC COEFFICIENTS 
C NOW INCLUDES ZERO FORWARD SPEED COEFFICIENTS 


IMPLICIT NONE 

INTEGER? JJJ,II,J 

INTEGER*2 IDTAR 

REAL*8 TEMP1,TEMP2,TEMP3, TEMP 

REAL*8 UPFRAC,LOFRAC,URACER,URACERO,URACES,URACESO 
REAL*8 DRE3,DSPE3,DSSE3,DRE,DSPE,DSSE 

REAL*8 DRHS1,DRHS2,DRHS3,DRHS4,DRHS5,DRHS6 

REAL*8 ZRHS1,ZRHS2,ZRHS3,ZRHS4,ZRHS5,ZRHS6 

REAL*8 BRHS1,BRHS2,BRHS3,BRHS4,BRHS5,BRHS6 


REAL*8 DD(1800) 
COMMON /COMDD/ DD 


REAL*8 URUR,USUS,DELRCR,DELRCS,ACOF(8),BCOF(8),acp(4),bcp(4) 
REAL*8 DSP,DSS,DR,DB,UC,DT,RHS(6),E(6),P(6),QPROP, TA, YG, YB 
REAL*8 PHISTERN, VATS, WATS,ETA,CETA,ETAMI,BETAS 
REAL*8 KBETA,WT,BOUY 

REAL*8 ARAPPHI2,ARAPPHI4, CROSSLIM, SINTHETA 

REAL*8 COSTHETA, SINPHI, COSPHI 

REAL*8 UU, VV,WW,PP,QQ,RR,PR,PQ,QR,UV,UW,UP 

REAL*8 UQ,UR, VATTAIL,UUDB 

REAL*8 VP,VR,WP,WQ,SORVVWW,VSQRVVWW 

REAL*8 WSQRVVWW,CROSS(6), VORT(6) 

REAL*8 DRDR,DSDS,UUDS,UUDR,UABW,UABV,WABW 

REAL*8 QABQ, VABV,RABR,FXP 

REAL*8 PABP,DPDP,UUDP,UT,BETA,ALPHA,WR 

REAL*8 ARRAY(6), WTINH2O,DBDB 


REAL*8 XUDOT,XQQ,XRP,XVR,XWQ,XVV,XDRDR,XDSSDSS 

REAL*8 XDSPDSP,XDBDB 

REAL*8 XW.XQ,XWDSP,XWDSS,XODSS,XQDSP,XVDR,XRDR,XWW,XRR 
REAL*8 XDRE2,XDSSE2,XDSPE2 

REAL*8 YVDOT, YRDOT, YRU, YVU, YVABVU,YVV 

REAL*8 YDRYDRN YRABRYVABV 

REAL*8 YPU, YWP, YPABP,YPQ, YSTAR, YVABVR, YHVABV, YDRE, YDRE3 
REAL*8 ZWDOT,ZQDOT,ZQU,ZWU,ZABWU,ZWW 

REAL*8 ZDSS,ZDSSN,ZOABQ,ZWABW 

REAL*8 ZVP,ZSTAR,ZPR,ZDSP,ZDSPN,ZHWABW,ZDB 

REAL*8 ZDSSE,ZDSPE,ZDSPE3,ZDSSE3 
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REAL*8 KPDOT,KVDOT,KRDOT,KQR,KPU,KRU,KWP,KVU,K4S 
REAL*8 K8S,KLKDSP,KDSS 

REAL*8 KPABP,KWR,KDR,KDRN,KSTAR 

REAL*8 KDRE,KDRE3,KDSSE,KDSPE,KDSSE3,KDSPE3 

REAL*8 MQDOT,MWDOT,MQU,MWU,MWABWR.MABWU,MWW,MDSS 
REAL*8 MDSSN,MQABQ,MWABW 

REAL*8 MRP,MSTAR,MVP,MDSP,MDSPN,MHWABW,MDB 
REAL*8 MDSSE,MDSPE,MDSSE3,MDSPE3 

REAL*8 NRDOT,NVDOT,NRU,NVU,NVABVR,NABVU,NVV,NDR 
REAL*8 NDRN,NVABV,NRABR 

REAL*8 NWP,NPQ,NSTAR,NPU,NHVABV,NDRE,NDRE3 

REAL*8 IYY,IZZ IXX,SLEN,MASS,ZB,ZG,XB,XG,XSTERN 

REAL*8 XFORWARD,KTAINF,UAINF,IXZ,IXY,IYZ, YPDOT,NPDOT 
EQUIVALENCE (DD( 1), DT)(DD(10) TA) 

EQUIVALENCE (DD( 11), JJJ), (DD( 16), BETA) 

EQUIVALENCE (DD( 15), ALPHA) 

EQUIVALENCE (DD( 17), FXP ),(DD( 18), QPROP) 
EQUIVALENCE (DD( 19), ETA ),(DD( 20), СЕТА) 
EQUIVALENCE (DD( 55), E(1)),(DD( 61), Р(1)) 

EQUIVALENCE (DD( 67), DSS),(DD( 68), DSP) 
EQUIVALENCE (DD( 69), | DR),(DD( 71), RHS(1)) 
EQUIVALENCE (DD(179), DB) 

EQUIVALENCE (DD(260),CROSS(1)),(DD(266), VORT(1)) 
EQUIVALENCE (DD(272),VATTAIL) 


EQUIVALENCE (DD(301),XUDOT )(DD(303,XW  ) 
EQUIVALENCE (DD(305,XQ  )(DD(307,XQQ ) 
EQUIVALENCE (DD(309,XRR  )(DDGlII)XRP ) 
EQUIVALENCE (DD(313),XVR  J(DD(315,XWQ  ) 
EQUIVALENCE (DDG17),XVV  )(DD(319,XWW  ) 
EQUIVALENCE (DD(321),XDBDB )(DD(323,XDRDR ) 
EQUIVALENCE (DD(325),XDSSDSS )(DD(327),XDSPDSP ) 
EQUIVALENCE (DD(329),XDRE2 )(DD(331),XDSSE2 ) 
EQUIVALENCE (DD(333),XDSPE2 ) 

EQUIVALENCE (DD(1703),XWDSP )(DD(1705),XWDSS ) 
EQUIVALENCE (DD(1707),XQDSP_),(DD(1709),XQDSS ) 
EQUIVALENCE (DD(1711),XVDR_),(DD(1713),XRDR_ ) 


EQUIVALENCE (DD(338),YVDOT_ ),(DD(340),YPDOT ) 
EQUIVALENCE (DD(342),YRDOT _),(DD(344), YPABP__ ) 
EQUIVALENCE (DD(346),YPQ __ ),(DD(348),YRU _ ) 
EQUIVALENCE (DD(350),YPU _ ),(DD(352),YVU _ ) 
EQUIVALENCE (DD(354), YVABVR_ ),(DD(356), YVABVU ) 
EQUIVALENCE (DD(358),YVV __ ),(DD(360),YDR __ ) 
EQUIVALENCE (DD(362),YDRN __),(DD(364), YRABR_ ) 
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EQUIVALENCE (DD(366) YVABV )(DD(368,YWP  ) 
EQUIVALENCE (DD(370), YSTAR ),(DD(372), YHVABV ) 
EQUIVALENCE (DD(374) YDRE3 )(DD(376),YDRE ) 


EQUIVALENCE (DD(390),ZWDOT )(DD(392),ZQDOT ) 
EQUIVALENCE (DD(394),ZQU _ ),(DD(396),ZVP _ ) 
EQUIVALENCE (DD(398,ZSTAR )(DD(400,ZWU  ) 
EQUIVALENCE (DD(402),ZABWU )(DD(404,ZWW ) 
EQUIVALENCE (DD(406),ZDB ),(00(408),70$$ ) 
EQUIVALENCE (DD(410,ZDSSN )(DD(412,ZDSP ) 
EQUIVALENCE (DD(414,ZDSPN )(DD(416),ZQABQ ) 
EQUIVALENCE (DD(418),ZWABW ),(DD(420)ZHWABW ) 
EQUIVALENCE (DD( 422),ZDSSE3 ),(DD( 424),ZDSPE3 ) 
EQUIVALENCE (DD( 426),ZDSSE_),(DD( 428),ZDSPE ) 
EQUIVALENCE (DD(1717),ZPR _ ) 


EQUIVALENCE (DD(440),KVDOT_ ),(DD(442),KPDOT_ ) 
EQUIVALENCE (DD(444),KRDOT _ ),(DD(446),KQR _ ) 
EQUIVALENCE (DD(448),KPABP_ ),(DD(450),KPU _ ) 
EQUIVALENCE (DD(452),KRU __),(DD(454),KWP__) 
EQUIVALENCE (DD(456),KSTAR__),(DD(458),KVU _ ) 
EQUIVALENCE (DD(460),KI ),(20(462),К4$ ) 
EQUIVALENCE (DD(464),K8S  )(DD(466,KDR  ) 
EQUIVALENCE (DD(468),KDRN ),(DD(470,KDSS ) 
EQUIVALENCE (DD(472),KDSP ),(00(474)КОВЕ ) 
EQUIVALENCE (DD(476)KDRE3 )(DD(478)KDSSE ) 
EQUIVALENCE (DD(480),KDSPE )(DD(482,KDSSE3 ) 
EQUIVALENCE (DD(484),KDSPE3 ) 

EQUIVALENCE (DD(1721,KWR ) 


EQUIVALENCE (DD(494,MWDOT ),(DD(496),MQDOT ) 
EQUIVALENCE (DD(498,MRP  )(DD(500,MQU ) 
EQUIVALENCE (DD(502),MWU __ ),(DD(504),MSTAR_ ) 
EQUIVALENCE (DD(506),MWABWR ),(DD(508),MABWU_ ) 
EQUIVALENCE (DD(510)MWW  )(DD(512)MDSS ) 
EQUIVALENCE (DD(514)MDSSN )(DD(516)MDSP ) 
EQUIVALENCE (DD(518),MDSPN ),(DD(520,MDB ) 
EQUIVALENCE (DD(522)MQABQ )(DD(524)MWABW ) 
EQUIVALENCE (DD( 526), MHWABW ),(DD( 528),MDSSE3 ) 
EQUIVALENCE (DD( 530),MDSPE3 ),(DD( 532),MDSSE ) 
EQUIVALENCE (DD( 534),MDSPE ) 

EQUIVALENCE (DD(1723),MVP _ ) 


EQUIVALENCE (DD(546),NVDOT ),(DD(548),NPDOT ) 
EQUIVALENCE (DD(550,NRDOT )(DD(552),NPQ ) 
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EQUIVALENCE (DD(554)NPU  )(DD(556)NRU ) 
EQUIVALENCE (DD(558)NSTAR )(DD(560)NVU ) 
EQUIVALENCE (DD(562)NVABVR )(DD(564)NABVU ) 
EQUIVALENCE (DD(566),NVV ),(00(568), МОК ) 
EQUIVALENCE (DD(570,NDRN ),(00(572),МУАВУ ) 
EQUIVALENCE (DD(574,NRABR )(DD(576),NWP ) 
EQUIVALENCE (DD(578),NHVABV ),(DD(580),NDRE3_ ) 
EQUIVALENCE (DD(582),NDRE_ ) 


EQUIVALENCE (DD(593),SLEN _ ),(DD(596),IYY ) 
EQUIVALENCE (DD(598),IZZ —)(DD(600),IXX ) 
EQUIVALENCE (DD(602),IXZ )(DD(604)IYZ ) 
EQUIVALENCE (DD(606),IXY —),(DD(608),MASS ) 
EQUIVALENCE (DD(610),WT __ ),(DD(612),BOUY ) 
EQUIVALENCE (DD(614),ZB ),(DD(616),2G ) 
EQUIVALENCE (DD(618),XB ),(DD(620),XG ) 
EQUIVALENCE (DD(622),XSTERN ),(DD(624),XFORWARD) 
EQUIVALENCE (DD(626),YG ),(DD(628),YB ) 
EQUIVALENCE (DD(629),WTINH2O ) 

EQUIVALENCE (DD(630),CROSSLIM) 

EQUIVALENCE (DD(631),LOFRAC. ),(DD(632),UPFRAC ) 


EQUIVALENCE (DD(848),KTAINF ),(00(849),0АІМЕ ) 
EQUIVALENCE (DD(850), URACERO ),(DD(851), URACESO ) 
EQUIVALENCE (DD(858),URACER ),(DD(859),URACES ) 
EQUIVALENCE (DD(866)DRE  )(DD(867),DSSE ) 
EQUIVALENCE (DD(868),DSPE ) 

EQUIVALENCE (DD(1758),ARRAY(1)),(DD(1764),ACOF(1)) 
EQUIVALENCE (DD(1800),botdep), (DD(1772),BCOF(1)) 
EQUIVALENCE (DD(1780),acp(1)), (DD(1788), bep(1)) 


C FOR EFFICIENCY THE FOLLOWING VARIABLES ARE CALCULATED 


С 


UU = E(1)*E(1) 
VV = E(2)*E(2) 
WW = E(3)*E(3) 
ОТ = DSQRT(UU+VV+WW) 
PP = E(4)*E(4) 
QQ = E(5)*E(5) 
RR = E(6)*E(6) 
PR = E(4)*E(6) 
PQ = E(4)*E(5) 
QR = E(5)*E(6) 
UV = E(1)*E(2) 
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UW - E(1)*E(3) 

UP - E(1)*E(4) 

UQ - E(1)*E(5) 

UR - E(1)*E(6) 

VP - E(2)*E(4) 

VR - E(2)*E(6) 

WP - E(3)*E(4) 

WQ = E(3)*E(5) 

WR = E(3)*E(6) 

URUR = URACER*URACER 
USUS = URACES*URACES 
DELRCR = URACER - URACERO 
DELRCS = URACES - URACESO 
VATS = E(2) + XSTERN*E(6) 
WATS = E(3) - XSTERN*E(5) 


IF( DABS(URACER) .GT. 0.1) THEN 

DRE = DR - DATAN2( VATS,URACER ) 
ELSE 

DRE = DR 
ENDIF 


IF( DABS(URACER) .GT. 0.1) THEN 
DSSE = DSS - DATAN2( WATS,URACES ) 
DSPE = DSP - DATAN2( WATS,URACES ) 
ELSE 
DSSE = DSS 
DSPE = DSP 

ENDIF 


DRE3 = DRE*DRE*DRE 
DSPE3 = DSPE*DSPE*DSPE 
DSSE3 = DSSE*DSSE*DSSE 
DRDR = DR*DR 

DBDB = DB*DB 

DSDS = DSS*DSS 

DPDP = DSP*DSP 

UUDR = UU*DR 

UUDB = UU*DB 

UUDP = UU*DSP 

UUDS = UU*DSS 

UABW = Е(1)*АВЅ(Е(3)) 
UABV = E(1)*ABS(E(2)) 
WABW - E(3)*ABS(E(3)) 
QABQ - E(5)*ABS(E(5)) 
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С 


VABV = E(2)*ABS(E(2)) 

PABP = E(4)*ABS(E(4)) 

RABR = E(6)*ABS(E(6)) 
SQRVVWW = DSQRT(VV+WW) 
VSQRVVWW = E(2) * SORVVWW 
WSQRVVWW -E(3)* SORVVWW 
WTINH20 =WT - BOUY 

SINPHI = DSIN(P(4)) 

COSPHI = DCOS(P(4)) 
SINTHETA = DSIN(P(5)) 
COSTHETA = DCOS(P(5)) 


IF( (DABS(VATS) .GT. 0.001) .OR. (DABS(WATS) .GT. 0.001)) THEN 
TEMPI - VATS*VATS 
TEMP2 - WATS*WATS 
TEMP3 = TEMP1+TEMP2 
ARAPPHI2 = TEMP3*2.*VATS*WATS/TEMP3 
ARAPPHI4 = TEMP3*4.*VATS*WATS*(TEMP2 - TEMP1)/(TEMP3*TEMP3) 
BETAS =DATAN2(DSQRT(TEMP3),E(1)) 
PHISTERN = DATAN2(WATS,-VATS) 
ELSE 
ARAPPHI2 = 0.0 
ARAPPHI4 = 0.0 


BETAS =0.0 
PHISTERN = 0.0 
ENDIF 


IF((DABS(E(3)) .GT. 0.1D-10) .AND. (DABS(E(1)) .GT. 0.01)) THEN 
ALPHA = DATAN2(E(3),E(1)) 

ELSE 
ALPHA = 0.0 

ENDIF 


IF( (DABS(E(2)) .GT. 0.1D-20) AND. (DABS(UT) .GT. .01)) THEN 
BETA = DASIN(-E(2)/UT) 

ELSE 
BETA = 0.0 

ENDIF 


C FOR EFFICIENCY CALL THESE SUBROUTINES ONLY ONCE 


С 


IF ( JJJ .EQ. 1) THEN 
CALL VORTEX 
CALL CROSFLOW 
CALL PROP 
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С 
DO 345 IDTAR = 1,10 
CALL TOWDARRY(IDTAR) 
345 CONTINUE 
ENDIF 
2 
СЕТА 
ЕТАМІ = ЕТА - 1.0 
C 
C CALCULATE RIGHT HAND SIDES OF EQUATIONS OF MOTION 
C 
C AXIAL FORCES ZERO SPEED MODEL 


C 
ZRHSI - 0.0 
С 
C LATERAL FORCE ZERO SPEED MODEL 
1 ZRHS2= YHVABV*VATS*DABS(VATS) 
c NORMAL FORCE ZERO SPEED MODEL 
1 ZRHS3 - ZHWABW*WATS*DABS(WATS) 
: ROLL MOMENT ZERO SPEED MODEL 
7 ZRHS4 = KPABP*PABP 
: PITCH MOMENT ZERO SPEED MODEL 
i ZRHS5 - MHWABW*WATS*DABS(WATS) 
E YAW MOMENT ZERO SPEED MODEL 
: ZRHS6 = NHVABV*VATS*DABS(VATS) 
C 


C CALCULATE DTRC HYDRODYNAMIC FORCES 
C HYDRODYNAMIC AXIAL FORCES 
С 
DRHS1= XQQ*QQ + XRR*RR + XRP*PR + XVR*VR + XWQ*WQ + XVV*VV 
1+ XW*UW + XWW*WW + XQ*UQ 
1 + (XDSSDSS*DSDS + XDSPDSP*DPDP + XDRDR*DRDR + 
1 + XDBDB*DBDB)*UU 
1 + (XWDSP*DSP + XWDSS*DSS)*UW + (XQDSP*DSP + XQDSS*DSS)*UQ 
| + (XVDR*UV + XRDR*UR)*DR 
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C HYDRODYNAMIC LATERAL FORCE 
C 
DRHS2 = YPABP*PABP + YPQ*PQ 
+ YRU*UR + YPU*UP + YVU*UV + YVABVU*UABV + YVV*VV + 
+ YWP*WP 
- (YDR + YDRN*ETAMI)*UUDR + YVABV*VABV + YRABR*RABR + 
+ YSTAR*UU 
+ YVABVR*VSQRVVWW 


мәй мә ээд і 


С 
C HYDRODYNAMIC NORMAL FORCE 
C 
DRHS3 = ZQU*UQ + ZWU*UW + ZABWU*UABW + ZVP*VP + ZSTAR*UU 
1 + (ZDSS + ZDSSN*ETAMI) * UUDS + (ZDSP + ZDSPN*ETAM1) * UUDP 
1 + ZDB*UUDB 
1 + ZQABQ*QABQ + ZWABW*WABW + ZPR*PR + 
1 + ZWW*DABS(E(3))*SQRVVWW 
C 
C HYDRODYNAMIC ROLL MOMENT 
C 
KBETA = BETAS*BETAS*(UU + VATS*VATS + WATS*WATS) 
DRHS4 = KPU*UP + KRU*UR + KQR*QR + KWP*WP = 
1 + ((KDR + KDRN*ETAM1)*DR + KDSS*DSS + KDSP*DSP + KSTAR)*UU 
1 + KWR*WR + KVU*UV + KPABP*PABP + KI*VATTAIL + ARAPPHI4*K4S 
С 1+KBETA*(K4S*DSIN(4.*PHISTERN) + K8S*DSIN(8.*PHISTERN)) 
С 
C HYDRODYNAMIC PITCH MOMENT 
C 
DRHS5 = MQU*UQ + MQABQ*QABQ + MWABW*WABW 
1 + MVP*VP + MRP*PR + MSTAR*UU + MDB*UUDB 
1 + MWU*UW + MWABWR*WSQRVVWW + MABWU*UABW + 
1 + MWW*DABS(E(3))*SQRVVWW 
1 + (MDSS + MDSSN*ETAM1)*UUDS + (MDSP + MDSPN*ETAM1)*UUDP 
С 
C HYDRODYNAMIC YAW MOMENT 
C 
DRHS6 = NPQ*PQ + NPU*UP + NWP*WP + NSTAR*UU 
| + NRU*UR + NVU*UV + NVABVR*VSQRVVWW + NABVU*UABV + 
1+ NVV*VV 
1 + (NDR + NDRN*ETAMI)*UUDR + NVABV*VABV + NRABR*RABR 
С 
C CALCULATE BODY FORCES, AND SUM FORCES INDEPENDENT OF ALPHA 
OR BETA 
C USED FOR BOTH LOW AND HIGH ANGLE OF ATTACK MODELS 
C 
C BODY AXIAL FORCE 
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С 


BRHSI - -MASS*( WQ - VR - XG*(QQ + RR) + ZG*PR) 

| - WTINH2O*SINTHETA + ЕХР + ARRAY(1) 

] + XDRE2*DRE*DRE*URUR 

1 + (XDSSE2*DSSE*DSSE + XDSPE2*DSPE*DSPE)*USUS 


C BODY LATERAL FORCE 


C 


© 
С 


BRHS2 = -MASS*(UR - WP + XG*PQ + ZG*QR ) 

| + WTINH20*COSTHETA*SINPHI 

1 + CROSS(2) + VORT(2) + ARRAY(2) 

1 + YDRE*( DELRCR*(-XSTERN*E(6) - E(2)) + DR*URUR) 
1 + YDRE3*DRE3*URUR 

1 + YDRE3*DRE3*URUR 


C BODY NORMAL FORCE 


C 


C 


BRHS3 = -MASS*(VP - UQ + XG*PR - ZG*(PP + QQ)) 
I + WTINH2O*COSTHETA*COSPHI 

1 + CROSS(3) + VORT(3) + ARRAY(3) 

1 + ZDSSE*( DELRCS*(XSTERN*E(5) - E(3)) + DSS*USUS) 
1 + ZDSPE*( DELRCS*(XSTERN*E(5) - E(3)) + DSP*USUS) 
1 + (ZDSPE3*DSPE3 + ZDSSE3*DSPE3)*USUS 


C BODY ROLL MOMENT 


C 


C 


BRHS4 = MASS*ZG*(UR - WP) 

1 - (ZG*WT - ZB*BOUY)*COSTHETA*SINPHI 

I + (YG*WT - YB*BOUY)*COSTHETA*COSPHI 

1 + QPROP + VORT(4) + ARRAY(4) 

1 + KDRE*( DELRCR*(XSTERN*E(6) - E(2)) + DR*URUR) 

| + KDRE3*DRE3*URUR 

1 + KDSSE * (DELRCS*(XSTERN*E(5) - E(3)) + DSS*USUS) 
1 + KDSPE * (DELRCS*(XSTERN*E(5) - E(3)) + DSP*USUS) 
I + (KDSPE3*DSPE3 + KDSSE3*DSPE3)*USUS 


C BODY PITCH MOMENT 


С 


BRHSS = (-IXX + IZZ)*PR + MASS*(-ZG*(WQ - VR) + XG*(VP - UQ)) 
1 - (XG*WT - XB*BOUY)*COSTHETA*COSPHI 

1 - (ZG*WT - ZB*BOUY)*SINTHETA 

1 + CROSS(5) + VORT(5) + ARRAY(5) 

I + MDSSE * (DELRCS*(XSTERN*E(5) - E(3)) + DSS*USUS) 

1 + MDSPE * (DELRCS*(XSTERN*E(5) - E(3)) + DSP*USUS) 

I + (MDSPE3*DSPE3 + MDSSE3*DSPE3)*USUS 
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C BODY YAW MOMENT 


e 


C 


BRHS6 - (-IYY + IXX)*PQ - MASS*(XG*(UR - WP)) 

1 + (XG*WT - XB*BOUY)*COSTHETA*SINPHI 

1 + (YG*WT - YB*BOUY)*SINTHETA 

1 + (XG*WT - XB*(BOUY))*COSTHETA*SINPHI 

I + (YG*WT - YB*(BOUY))*SINTHETA 

1 + CROSS(6) + VORT(6) + ARRAY(6) 

1 + NDRE*(DELRCR*(-X STERN*E(6) - E(2)) + DR*URUR) 
1 + NDRE3*DRE3*URUR 


C SUM ALL FORCES 


5 


RHS(1) = BRHS1 + LOFRAC*DRHS1 + UPFRAC*ZRHS I 
RHS(2) = BRHS2 + LOFRAC*DRHS2 + UPFRAC*ZRHS2 
RHS(3) = BRHS3 + LOFRAC*DRHS3 + UPFRAC*ZRHS3 
RHS(4) = BRHS4 + LOFRAC*DRHS4 + UPFRAC*ZRHS4 
RHS(5) = BRHS5 + LOFRAC*DRHS5 + UPFRAC*ZRHS5 
RHS(6) - BRHS6 + LOFRAC*DRHS6 + UPFRAC*ZRHS6 


RETURN 
END 


87 





Appendix B. Trajectory input files. 


This appendix presents the geometry limitation file and trajectory command files used to 
generate the simulation trajectories. 
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Geometry Limitation File GEOM.LIM 


93.57 16.50 0. ХС YC ZC [FT] 

85.275 15.04 0. XMIN YMIN ZMIN [FT] 

93.576 16.50 0. XMAX YMAX ZMAX [FT] 

8.00 16.5 0.8125 10.5 D2 HRAD VRAD SLEN [FT] 

0.1 1.625 2.0 1.47 VHOST [FT/SEC] DIA [FT] RHOC [SLUG/FT^3] 
ATLMAX 
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Traj 1 


2.0 0.04 VELOCITY AND ACCELERATION EXPONENTS 
XVEL,YVEL,ZVEL CORRECTION COMMANDS [FT/SEC] [SEC] 
PIA PIT P2A P2T РЗА РТ 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


ROLL,PITCH, YAW RATE COMMANDS [DEG/SEC] [SEC] 


P4A 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
20.00 
30.00 
35.00 
40.00 

799.00 
899.00 
900.00 
990.00 
999.00 


0.00 0.00 0.00 0.00 


0.00 600.00 
0.00 700.00 
0.00 750.00 
0.00 800.00 
0.00 850.00 
0.00 860.00 
0.00 870.00 
0.00 880.00 
0.00 900.00 


0.00 600.00 
0.00 700.00 
0.00 720.00 
0.00 750.00 
0.00 800.00 
0.00 810.00 
0.00 820.00 
0.00 850.00 
0.00 900.00 


РАТ PSA PST РбА  P6T 
0.00 0.00 0.00 0.00 0.00 


600.00 
700.00 
750.00 
800.00 
900.00 
999.00 
999.00 
999.00 
999.00 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


600.00 
700.00 
750.00 
800.00 
850.00 
860.00 
870.00 
880.00 
900.00 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


20.00 
25.00 
28.00 
750.00 
800.00 
810.00 
820.00 
850.00 
900.00 
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Тга) 2 


20 004 VELOCITY AND ACCELERATION EXPONENTS 
XVEL,YVEL,ZVEL CORRECTION COMMANDS [FT/SEC] [SEC] 


РТА 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


ROLL,PITCH, YAW RATE COMMANDS [DEG/SEC] [SEC] 


РАА 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


PIT 
0.00 
20.00 
30.00 
35.00 
40.00 
799.00 
899.00 
900.00 
990.00 
999.00 


P2A PE 


РЗА BST 


0.00 0.00 0.00 0.00 


0.00 600.00 
0.00 700.00 
0.00 750.00 
0.00 800.00 
0.00 850.00 
0.00 860.00 
0.00 870.00 
0.00 880.00 
0.00 900.00 


0.00 600.00 
0.00 700.00 
0.00 720.00 
0.00 750.00 
0.00 800.00 
0.00 810.00 
0.00 820.00 
0.00 850.00 
0.00 900.00 


РАТ PSA PST РбА P6T 
0.00 0.00 0.00 0.00 0.00 


600.00 
700.00 
750.00 
$00.00 
900.00 
999.00 
999.00 
999.00 
999.00 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


600.00 
700.00 
750.00 
$00.00 
850.00 
860.00 
870.00 
880.00 
900.00 


0.00 
-2.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


20.00 
25.00 
30.00 
35.00 
36.00 
$10.00 
820.00 
850.00 
900.00 
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Тга) 3 


2.0 0.04 VELOCITY AND ACCELERATION EXPONENTS 
XVEL, YVEL,ZVEL CORRECTION COMMANDS [FT/SEC] [SEC] 
PA PIT PA PT РЗА РЗТ 


0.00 
0.00 
0.00 
0.00 
-4.00 
-4.00 
-4.00 
-4.00 
-4.00 
-4.00 


ROLL,PITCH,YAW RATE COMMANDS [DEG/SEC] [SEC] 


РАА 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
10.00 
15.00 
20.00 
30.00 

799.00 
899.00 
900.00 
990.00 
999.00 


P4T 
0.00 
600.00 
700.00 
750.00 
800.00 
900.00 
999.00 
999.00 
999.00 
999.00 


0.00 
0.00 
0.00 


0.00 600.00 
0.00 700.00 
0.00 750.00 
0.00 800.00 
0.00 850.00 
0.00 860.00 
0.00 870.00 
0.00 880.00 
0.00 900.00 


PSA PST 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


600.00 
700.00 
750.00 
800.00 
850.00 
$60.00 
870.00 
880.00 
900.00 


0.00 


0.00 0.00 0.00 0.00 


0.00 600.00 
0.00 700.00 
0.00 720.00 
0.00 750.00 
0.00 800.00 
0.00 810.00 
0.00 820.00 
0.00 850.00 
0.00 900.00 


P6A P6T 


0.00 0.00 0.00 0.00 


0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 


20.00 
25.00 
28.00 
750.00 
800.00 
810.00 
820.00 
850.00 
900.00 


22 








Тга) 4 


Traj 4 uses Traj 1 with Vhost = 0.01 in Geom. lim 


93 








Appendix C. Model Fit Analysis of Variance. 
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Source Freedom Squares | 
| Modl | 1  .[0.04653676 | 0.046537 | 4.0022 | 0.0627 . 
| Emo | 16  [|o18605s5 | 0011622 | — 


MEL ”СЛСЭР? ЭЙ” ШҮРСРЭР A a 
| 34 | Мм | 1 | 0.13460215 | 0.134602 | 26603 | 0.1339 
| в | 16 050532 | 00506 | |  — 
тээг эт? лт й > НН о ЈА 
210 | Мода |  — 1 | 11940734 | 119407 | 257911 | «00001 | 
________-епи | 7 32408566 | 0490 | — 
| TRATA aao o A 
116 | Modi | 1 | 32888807 | 3.28888 | 889794 | <0.0001 | 
[ | mor | 16 | 05013961 | oo e aaa 
г засаа AAA SAA src 
PAR 


1.32 1 | 4.7664514 | 4.76645 | 166.9032 «0.0001 


0.9709780 0.02856 


НЭТ га USUS -— 
1.46 Model — 0.46515718 | 0.465157 31.1761 حال‎ 





0.14929 0.01492 


0.61436 ж E 


















| Model | 1 | 0.04653676 | 0.046537 | 4.0022 — 
[ Emo | 16 0186058 | oomes | — 
ош ao pezeser| | | 
PE MS E ы ла ЭР" 









==" 
[1 [1194054 | 19407 | 257911 | 500001 _ 
or | to [xen [0m | a 
бы Га” 1. 245495005 аас TE Об å] 
ма Г 1—— [338807 Зэн | нэн | ог 
[mor | 16 p esme | ome |__| 
— и e) | 
47664514 | 476045 | 1669032 | <1 | 
0, 9709780 | 0, 02856 
pos as sm J s s= 
| Мон | 1 1 0.46515718 | 0.465157 | 311761 | «00002 | 
[ 1o  o14900 | 00192 | | | 
mm QUEM E NP 


Table C-2 ANOVA Table for PTF effect on Launchway Moment 














95 


14 5iNS VIOLE 
45 “тн 1115 [ 
1/99 22527-200 wie 1 




















4. иены 
š ТОВ 4 ҚАҚ 


m i 


Ап 
PAPO EL IC ICH 
ИШКЕН зон 

4 NN GS) qa ye 
% TUM 4 ЫРЫ 


ДАНС 
PT 
Ths “НЭГ $c "Lax 
AAR PA KAN y» 
۹ LI Ұ 
Хад АА ГЕЗ 


ЕЗИ: 
Skre Чың 4 
4.8. 

РАЈ H 


Ма 77, 


међа: 
OA ONA) 
FS RA LLA AE 
igh ЕТСЯ 
ИСН 
үн, * =H; ot 
Fal, Ar 


шош 
T hør vi 

bi ek ope ib 
T %. 4% 1 
ПЕ 4 
ХАЛЗ: vea ser 
ЖЕТЕГІ D 
ТҮГҮЛ ЖӨ (12 
палатаны 


eti 
ч к. д 


КЮ 


v 
Py ksta 
. 


Дэ IU 


Dr НЭЭ? 


РТ ТУС 


"M 

55 4 vike ғ. 

илгээ saat ares] one? 

С z ha Ч 

Ч ¿AUN ТЭЭЛ 

84646467 НЕТ 

ОДА ОР МАС И 

A ol » 

ti na” rå 
8 КГ: rt 


AA MT ғ på 
^4۹ 


Т 


Я 


логе 
КОЛЕ РЕКЕ 
дээ 
... 


.. ер З Dor 
DL LET 
* a farata 


[NE 


Vs d yx 


119 | 
МИ? 
21 F 


Ж УА 

A 

TS E 
e HERE ay 


БИ м 
"t wf AI 


P 


M AM. 


«16 ТУХ) 
ҚАЛЫ ыы 


ur 


[LIAE 
m PEN Ma! Қ 
b pi (77 


му АТ 
МА 


"ap Cala 2172 
Ж” Чи 


114 их 
ae А А 
АЛ 
ona М; (22 Å 

"У 


ИЛ 


iets 


*, 
d PALS : 
Я 
vente 
«lie 


ДОС РРА ЕН 
АДЕ 


"пи OE M ACO Wo 
AJENO bea gen i у M av 
- 


КДУ i n 


СЕГТОСО”Ж 


SM 
må Гын, a. 1 
e 


با 
Tan АНЕ‏ 


УГ! moe v. 
DU Ч: 


ЕРА 

Ши 

Iu 

tie EE 

Li ` 

^e - NI > 
EU 
|) ' > 
B uy кв 

од | 

ir Ч reg eta 


ШЕ 


nm vp ty ©, D 


N 7 TELAT TL] 


Ls ç r h ЖИ oe Зө 


9^4 * 3 «. PT 


џи "i m 


KERLE НАМ ПИР 


LELEL TIELA Ty 


LETTA SA 
n Do VL Le TS ҮТ? 
© "ir ГАЛЛЕ УЫ . 
i ا‎ ty НО 1. 438, 
PES Бли А ДОМЕ ДОМУ 
1 t 279 
ITI TR T 
Ы UO 
D 4 алата | | veh) 
% НОЈ) Pr adt git 3 
bh Geta aot m. 
TON [ Ir А tar AT TD мын 
N ГҮ .. им, D 
кес LP | 
1,6/5,4 мал а 
X Hi re“ 
М 


EA 


Дэг 


505172 


wi 
LU 2 ЖЫ 
НЭГ” 

ПОО У) 
Md And ҮШ [22 
D OPERI | 
PE m | 


б“. 


k. 
FT 


(1,2 


sy P 
"м1 " orm зе ns H 


AF XP 
%, M 
ЛЕЛЕ Op 
Jen DO ES А СИМ?” 
ЭГ б 
TE Em TEN 


A EP 
21113 HOT 
Г al Ls 
syte 4 på å ! 4 SA] 
Hon aat, z 2, PT 
4 М. 
se STE PG GAEL 
OE GE АСҚА 
Ци E 
ЦЭР 
4747.40 ..4 
LTE ID P 
(РРА 


Ын. 

ТИГЕТ 
ЦАЖ 121 
LUNES C] 
Lr zT 


ULL T М .% 
CET ORE r 
шинэ 
ЖҰТ 


DD. EET 


8404 
ч 


x 


.... 


rr ma 


4 
КІН; 


Tn 
ЯНГ 
T H 
. TE 
wA | 
¡Husa e 


DADA i 

Н АЛТ 

И tee k 
EEUU С 


H П 
СИЛУ 


Цэг: 
ЫНТЫ 
h жин HU 
^um HERE , 
“ 


TEES vio 
p" шил 


«унэ 
“ " 
Но ар Тір 
Ay LG 
ЦЭН ба 
Ч AP ED M ЛИЕ; 
4) Å у/о«& 12 e 
satu Н, “Ч К 
AJ Ud SLE TUR 
(44.3 XY H КҮЛГЕ R 
Ж мм 
ННН S 
ХАЯХ 8.4. 
TE “е Т 
ev. LS AU de bcd e 
Хэ! LU 
E LU 


. пи 
ye 112127 
i EHE ҚАН A 

Я? ЕЕ 
ТЕТ V 0 de 
УЧНИРЧЫН 


SØNN 
дор шал, 
PP ' p 
Em 


Г wie uv 
eras $ 
тд 
КД ХЕ 3 ERARI 1i Hs 
A) НИЙ ІЛЕ АНЫН 
ҒЫН ТА? 
ei wt Ad ru village 10 
ККД ОНИ 
А imu ^W HCM ad E Tor 
Ц С ИГЫ DOE 
у’ МИ ТН: ДЕА 
us ig KUTTER OI TE 4. К 
22 rn С Qe tat ин КУЕ P 
тийн 4 К ПАР И ИГ 
jn utt НҮД НЬ А НП TA 579 1⁄4, 49; 
Н ue 


t п КИЕ " 
9 дүн nhe ОЮ) ~ 
мать WAS: LO) 
13 


p 
мэ Бэ АЧ Ч 


№ f! 
СИЈА 
ШУ ni" 


БН a? at 
П 


Ц 
ёс D 
Set 
ЧЭ , n E 


i M un 


1 
4. a 
LATE 1 
ТН 1 0 


IG Н 


EE 


ЕТТІ DET 
11171 
KE 


ta. 


a q) 


nn 
ЦОО 
ray 
MLE d a 

7%, A en. ы-ы 
P [ ‘ 


IURE УЛ! 


ЖА 


1241 1 M, ЖАТЫ 
КИЛ ат РАЧЕ 


NUT БОГ 
маза ч 


АТЫ PL ат 
цг 
1 


(КИТЕ 


“ 9 
CELE ETET ETE" E 
en П 


СИИ. 
‘ 


аа 
t ditag t зан 


"р А 


Ганц 
Алымы 


LU] "n 


fe gf af u 
3 

å du TEE 

Je 


" и’ 
из 
М 


å Hi М 
nt 
Н {) ЭН 
kite ke 
JN TUE R 
ү (4 Ч еп СЫ и, 
НА АЛА И 
д лай үч 
112 ha: Á 
Eur 
Un à 
BM TM 
få шу oa 


Pu HER АВЧ 
DL 
h А D 
e 
И 216 ТИР 


малыг N^ 
Pat 


mo DUE 
СУ SA 
% n Л b Ex RS 


la 


Le PEE vs 


CLEF STEE Er) 


и 
ПАКИ ГЕ 
DEREN ПР 


САР 
ХЭН 


` 
IM 


45 у va 
u ON 
174 


OLE 
PD 


eit 
аа n 
А ОК 
LITI 


pa 
ог 
= ГЫС På 
ШИК ЕГ 
TT 


2112) 
БАЛТ 





